Geometrical Aspects in Optical Wavepacket Dynamics 
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We construct a semiclassical theory for propagation of an optical wavepacket in non-conducting 
media with periodic structures of dielectric permittivity and magnetic permeability, i.e., non- 
conducting photonic crystals. We employ a quantum-mechanical formalism in order to clarify its 
link to those of electronic systems. It involves the geometrical phase, i.e.. Berry phase, in a nat- 
ural way, and describes an interplay between orbital motion and the internal rotation. Based on 
the above theory, we discuss the geometrical aspects of the optical Hall effect. We also consider 
a reduction of the theory to a system without periodic structure and apply it to the transverse 
shift at an interface reflection/refraction. For generic incident beams with elliptic polarizations, an 
identical result for the transverse shift of each reflected/transmitted beam is given by the following 
different approaches; (i) analytic evaluation of wavepacket dynamics, (ii) total angular momentum 
(TAM) conservation for individual photons, and (iii) numerical simulation of wavepacket dynamics. 
It is consistent with a result by classical electrodynamics. This means that the TAM conservation 
for individual photons is already taken into account in wave optics, i.e, classical electrodynamics. 
Finally, we show an application of our theory to a two-dimensional photonic crystal, and propose 
an optimal design for the enhancement of the optical Hall effect in photonic crystals. 

PACS numbers: 03.65.Vf, 42.15.-i, 42.15.Eq, 42.70.Qs, 



I. INTRODUCTION 



The geometrical phase known as the Berry phase [1] 
has been attracting extensive interests in various fields, 
e.g., optics, molecular physics, nuclear physics, and con- 
densed matter physics [2, 3]. In particular, in con- 
densed matter physics, important roles of the geomet- 
rical phase in electronic transport phenomena have been 
intensively studied in the past several years, and great 
strides has been made both in theoretical and experi- 
mental researches. Although a hint of the Berry phase 
was recognized long time ago as the anomalous velocity 
in fcrromagnets which leads to the anomalous Hall ef- 
fect [4, 5], it was only after the discovery of the quantum 
Hall effect that the role of the Berry phase in electron 
transport began to be recognized. In the quantum Hall 
system under strong magnetic field, the Hall conductance 
was related to the topological integer, i.e., Chern num- 
ber [6-8]. A recent development is the finding that the 
Berry phase structure is a fundamental characteristic of 
the Blocli wavef unctions even in ordinary systems. From 
this viewpoint, the similarity between the anomalous and 
quantum Hall effects has been revealed [9-11]. The spin 
Hall effect based on the geometrical mechanism are also 
proposed recently and opened a new stage of spintron- 
ics [12, 13]. All of these effects are understood from the 
concept of the generalized anomalous velocity due to the 
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Berry phase. In other words, a trajectory of an electron 
is affected by the Berry phase. 

In optics, one can also find some phenomena related to 
the Berry phase. A change of polarization of light dur- 
ing propagation such as in helically wound optical fibers 
found in the early days [14-16] has been related with the 
Berry phase [17-19]. It's influence on the trajectory of 
light has been studied recently by deriving a set of semi- 
classical equations of motion [20]. From this viewpoint, 
there are several optical phenomena which are now inter- 
preted as a change of light trajectory due to the Berry 
phase. One is a transverse shift in reflection/refraction 
at an interface between two homogeneous media [21-28], 
(This effect in the case of internal total reflection is called 
as Imbert-Fedorov shift.) The other is a rotation of the 
beam inside an optical fiber, which is sometimes called as 
an optical Magnus effect [29-31]. These phenomena can 
be coined as the optical Hall effect, because of the sim- 
ilarity to the topological Hall effects [6-13] in electronic 
systems. 

Attribution of these optical phenomena to the Berry 
phase is not merely a re-interpretation, but also can open 
a new frontier for novel phenomena. The present authors 
[20] proposed that in photonic crystals the optical Hall 
effect is enhanced by an order of magnitude than the 
above-mentioned examples. It is inspired by its electronic 
counterpart; the topological Hall effects are known to be 
enhanced by periodic potentials, particularly when two 
bands come close in energy. Thus by designing a photonic 
crystal to have near-degenerate bands, the predicted shift 
of a light beam is large enough to be observable in exper- 
iments [20] . To calculate and design such photonic crys- 
tals in a quantitative way, a semiclassical Berry-phase 
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theory of optics in such photonic crystals is caUcd for. 
For this purpose, the approach by the variational princi- 
ple [32-35] , which we take in our previous [20] and present 
papers, works better than other approaches [30, 31]; it is 
because the approaches in Refs. [30, 31] use an eikonal ap- 
proximation, by which it is rather difficult to fully incor- 
porate the vectorial nature of the electromagnetic waves 
for generic cases like photonic crystals. 

In the previous work [20] , we have briefly reported the 
essence of the optical Hall effect and the mechanism of 
its enhancement in artificial crystals called photonic crys- 
tals, i.e., systems with periodic structures of dielectric 
permittivity /magnetic permeability [36]. In the present 
paper, we construct a semiclassical theory of an opti- 
cal wavcpackct (or a photon wavcpackct) in full detail 
by keeping its close connection to a theory of an elec- 
tron wavepacket. It incorporates the Berry phase in a 
natural way. The main focus of the present paper is to 
present basics to the extended geometrical optics appli- 
cable to photonic crystals. This class of artificial crystals 
are attracting great interests as new optical materials. 
Photonic crystals can be designed to have a desired band 
structure, with an aid of first-principle numerical calcula- 
tions, which enables the control of many novel properties 
of lights [36, 37] . To serve for such purposes, our theory 
is presented in a transparent way suitable for such nu- 
merical calculations. As is briefly presented in our pre- 
vious work [20], the effect of the geometrical phase on 
an optical wavepacket can be incorporated in the same 
manner as that in electronic systems. These generalized 
equations of motion correctly describe the interplay be- 
tween the orbital motion and the internal rotation, e.g., 
polarization, of wavepackets. The effect similar to the 
electrical Hall effects in condensed matter is derived in 
photonic systems with periodic structures. Indeed, our 
equations of motion are analogous to the semiclassical 
equations of motion for electron wavepackets in solids 
[34, 35]. However, the latter basically considered spin- 
less electrons, and in the case of optics, the polarization 
degrees of freedom has to be taken into account, where 
the Berry connection is non-Abelian in general. In this 
sense, an optical wavepacket is more similar to a spinful 
electron wavepacket. 

Below, for simplicity, we focus on a light propagating 
in a non-conducting medium in which there is neither 
electric nor magnetic order, i.e., the dielectric permittiv- 
ity and the magnetic permeability are symmetric tensors. 
Also their frequency dependences are neglected for sim- 
plicity. These conditions ensure the equation of continu- 
ity of electromagnetic energy [38] and we can construct 
the unitary theory of electromagnetic field. Based on 
this theory, the semiclassical equations of motion can be 
derived on an equal footing with electronic systems in 
which the semiclassical equations of motion arc derived 
from quantum mechanics. In order to stress the analogy 
between electronic and photonic systems, we formulate 
a theory for Bloch states of electromagnetic field in a 
quantum-mechanical formalism. Although we focus on 



the unitary theory in this paper, its extension to a non- 
unitary version for systems with electric/magnetic order 
and conducting systems would give some insights to the 
interesting phenomena and proposals, e.g. the photonic 
Hall effect in a scattering media subject to an external 
magnetic field [39, 40], the magnetically induced deflec- 
tion due to the Pitaefskii magnetization [41-46], the one- 
way waveguide of edge states in magnetic photonic crys- 
tals [47, 48], and Lorentz force on the light due to the 
toroidal moment [49]. 

The reduction to a system without periodic structure 
is straightforward. Indeed, in the previous work [20], we 
have presented a simple application of our theory to the 
transverse shift in the reflection/refraction at an inter- 
face, and found that this shift is governed by the conser- 
vation of total angular momentum (TAM) for reflected 
and refracted photons individually. We have also nu- 
merically demonstrated the validity of our theory for the 
case of an incident wavepacket with circular polarization. 
However very recently, the transverse shift evaluated by 
the conservation of TAM has been questioned [50] for the 
cases of incident wavepackets with elliptic polarizations. 
In this paper, we present an additional way of estimat- 
ing the transverse shift from the asymptotic form of the 
wavepacket, and also the numerical calculations for the 
generic polarized states. As an important consequence 
from the study on this issue, we find that an identical 
result for the shift of each beam is given by (i) analytic 
evaluation of wavepacket dynamics, (ii) TAM conserva- 
tion for individual photons in Ref. [20] and (iii) numer- 
ically exact simulation of wavepacket dynamics. This 
agreement in different approaches supports the validity 
of the present theory claiming that the transverse shift 
is governed by the conservation of TAM for individual 
photons. These results are also consistent with the more 
conventional approach based on classical electrodynam- 
ics [26, 27]. In other words, the TAM conservation for 
individual photons is already taken into account in wave 
optics, i.e., classical electrodynamics. 

For a broad readership, we divide the main contents 
into two sections; Sec. II is devoted to formalisms, ex- 
plaining in full details the derivations of the theory and 
the resulting formulae, while in Sec. Ill we focus on two 
applications of the theory: the transverse shift in inter- 
face reflection/refraction, and the optical Hall effect in 
a two-dimensional photonic crystal. Readers who are 
mainly interested in the applications can skip Sec. II and 
jump to Sec. III. For this purpose we make Sec. Ill to be 
self-consistent. 

The plan of this paper is as follows. In Sec. HA, 
an electromagnetic field in a non-conducting medium is 
quantized in the Hamilton- Jacobi formalism by introduc- 
ing the Dirac bracket for the constrained system. Some 
quantum operators for physical observables arc also pre- 
sented. Eigen states in a periodic system are discussed 
in Sec. II B for the application to a photonic crystal. In 
Sec. II C, we consider a perturbed modulation superim- 
posed on a background periodic structure and discuss 
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corrections for the cigcn states and the expectation val- 
ues of physical observables for an optical wavepacket. 
The equations of motion are derived taking into account 
the Berry phase and the perturbed modulation. An ap- 
plication of our theory to reflection/refraction problem 
at a flat interface is discussed in Sec. Ill A by reducing 
the theory to the case without periodic structure. Re- 
cently some criticisms are raised against our approach to 
this reflection/refraction problem [50]. Remarks on these 
criticisms are presented in Sec. Ill B. In Sec. Ill C, we ap- 
ply our theory to a modulated two-dimensional photonic 
crystal and present some examples of Berry curvatures 
and internal rotations in a periodic system. 

Section IV is devoted to the discussion on the implica- 
tions of the present work to wider range of phenomena 
in physics. Related previous works are mentioned here. 



II. FORMALISMS 



Electromagnetic field in a non-conducting 
medium 



We consider an electromagnetic field in a non- 
conducting medium with a generic modulation but with- 
out electric nor magnetic orders, and begin with the fol- 
lowing Lagrangian, 

L =\j dr[E{r,t)-D{r,t)-H{r,t)-B{r,t)],{l) 



where 

D{r,t) = 7{r)E{r,t) = 7ir)[-dtA{r,t) -Vr(l>{r,t)], 

(2a) 

B{r,t) = *il{r)H(r,t) = x A{r,t). (2b) 

We take the unit in which h = c = 1 where c is the 

speed of light in vacuum. The medium where light 
propagates is treated as an insulating material, which 
is characterized by the dielectric permittivity e*(r) and 
the magnetic permeability Jl{r). These are assumed to 
be locally symmetric and real-valued tensors, and their 
frequency dependences are neglected. As mentioned in 
Sec. I, the equation of continuity for the electromagnetic 
energy holds under these conditions [38]. It should be 
noted that they arc the sufficient conditions but not the 
necessary conditions. The functional derivatives with re- 
spect to dtA(r,t) and dt(f>(r,t) determine the canonical 
momenta 7r(r,i) for A{r,t) and ^^{r^t) for (f>{r,t) re- 
spectively. The former gives the canonical definition of 
7r(r,f) as 7r(r,f) = —D{r,t), while the latter gives the 
constraint 7r0(r,t) = 0. Here wc introduce the Lagrange 
multiplier X^{r,t) for this constraint, and the Hamilto- 



nian is given by 



H = j drTT{r,t)-dtA{r,t)- L + J drX4r,t)Tr^{r,t) 
= Ho + j dr[-Tz{r)-Vr(f>{r)+X4,{r)TT4r)], (3a) 

Ho = ^ y dr[7r(r)e*~^(r)7r(r) 

+ [V, X A{r)])l-\r)[Vr. x A{r)]]. 



(3b) 



In order for the constraint tt^ ~ to be consistently 
satisfied, the following additional constraint is required. 



{TT4,{r),H}p 



-V,. • 7r(r) « 0, 



(4) 



where {• • • }p is the Poisson bracket. The symbol "w" 
means "weak equality" which is satisfied when all con- 
straints arc imposed [51]. When the Poisson brackets 
among a set of constraints and a Hamiltonian vanish on 
a constrained subspace, these constraints arc called first- 
class constraints by definition. On the other hand, when 
the Poisson brackets of these constraints among them- 
selves do not vanish even on the constrained subspace, we 
call them second-class constraints. In the present case, 
7r,^(r) and Vt- • 7r(r) commute with each other, and the 
commutation relation between V,, • 7v{r) and the Hamil- 
tonian generates no additional constraint. So the present 
system has two of first-class constraints. 



XI (r) = 7r0(r) « 0, 
X2{r) = Vr ■ 7r(r) w 0. 



(5a) 
(5b) 



In order to make a canonical formalism for such a con- 
strained system, all the first-class constraints arc trans- 
formed to be second class by introducing gauge fixing 
conditions. Here we take the following gauge conditions. 



X3(r) = (/)(r) « 0, 

X4(r) = Vr^r)A{r)^0. 



(6a) 
(6b) 



Then the commutation relations between the original 
constraints and the gauge conditions are represented by 



C{r,r') = {xa{r),X0{r')}p 

_ ( -C{ry) 



C(r,r') 



C(r, r') 
5{r - r') 







-Vre{r)Vr5{r - r') 



(7a) 



(•7b) 



Introducing Lagrange multipliers for the constraints 
including the gauge conditions, we redefine the Hamilto- 
nian as 



H = Ho + J drXir) ■ x{r), 



(8) 



where A(r) = [Xi{r), X2{r), Xs{r), X4{r)] and x(r) = 
[Xi{r),X2{r),X3{r),X4{r)]. The Lagrange multipliers 
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are determined by the conditions {x{r),H}p w and 
given by 

A(r) = - J dr'C-'{ry){x{r'),Ho}p, (9) 



where 



(lOb) 



g{r,r') 

and g{r,r') satisfies 

Vr^{r)Vr9{r,r') = Vr''?^{r')Vr'g{r,r') 

= -S{r-r'). (11) 

As a preparation for the quantum theory, we introduce 
the Dirac bracket defined by 

{F,G}j, 

= {F,G}p- J drdr'{F,x{r)}pC-\r,r'){x{r'),G}p. 

(12) 

Especially for A{r) and 7r(r), we obtain the following 

relation, 

{AXr),7r,(r')}D 

= S^.5{r - r') - ^ v;eJ,(r')Vj;,<?(r,r'), (13) 

k 

and this leads to the relation between the physical ob- 
servables, 

{B,(r),D,(r')}D = e,,feV^^(r - r'), (14) 

k 

where e^fe is the completely antisymmetric tensor defined 



using e 



xyz 



1. It is noted that this relation is the same as 
that in the vacuum, while that between E{r) and H{r) 
is not the case. The equations of motion are derived as 

dtA{r,t) = {A(r,t),i?}D« e-'(r)7r(r,0, (15a) 
dtTz{r,t) = {Tv{r,t),H}n 

= -Vr X [rt-\r)[Vr X A{r,t)]\ .(15b) 

The above equations are equivalent to the Maxwell equa- 
tions, 

dtDir, t) = X H{r, t), (16a) 
dtB{r, t) = -Vr X E{r, t), (16b) 
• D{r) = Vr ■ B{r) = 0. (16c) 

The present system is straightforwardly quantized by 
the identification as i{F,G}-D — > [P-iG]. Especially, the 
basic commutation relation is quantized as follows, 

[S,(r),D,.(r')] = i^e,,-,V^^(r-r'). (17) 



Here wc introduce some quantum operators which are 
useful to check the property of a wavcpackct. They are 
the Hamiltonian H, the center of the position TZ weighted 
by energy density, the energy current (the Poynting vec- 
tor) "P and the rotation of energy current JT, which are 
respectively defined by 



H = 


2j 


[ dr [E{r) ■ D{r) + H{r) ■ B{r)] , 


(18a) 


■R = 


h 


f drr [E{r) ■ D{r) + H{r) ■ B{r] 


1], (18b) 


V = 


2j 


f dr [E{r) x H{r) - H{r) x E{r 


)], (18c) 


J = 


o 1 


fdrrx [E{r) x H{r) - H{r) x 


E{r)]. 



(18d) 

It should be noted that the last two operators are differ- 
ent from the momentum and angular momentum opera- 
tors defined by 

P = dr [D{r) X B{r) - B{r) x D{r)] , (19a) 

J = ^jdrrx [D{r) x B{r) - B{r) x D{r)] , 

(19b) 

while V and S are conceptually close to P and J, re- 
spectively. This is because and S are not necessarily 
proportional to P and J. In other words, V and 
do not necessarily satisfy the algebra of the momentum 
and the angidar momentum. Therefore, in a system with 
translational and rotational symmetries, what should be 
conserved are P and J, rather than "P and Actually 
it has been experimentally confirmed that J is conserved 
in a dielectric medium with rotational symmetry [52] . In 
spite of these shortcomings of 7-' and JT", when a system 
has no continuous translational nor rotational symmetry, 
we focus on V and JT". This is because "P and JT" have 
relatively simple expressions even in a periodic system as 
shown in Appendix C. Especially, a part of J" suggests 
a close relation between the internal rotation and the 
Berry curvature in a photonic system, as well as in the 
quantum Hall system where the internal rotation of an 
spinless electron is originated by the cyclotron motion 
[34]. However, when a system has continuous transla- 
tional and rotational symmetries, we focus on P and J. 
This is the case in the reflection/refraction problem at a 
flat interface in Sec. HI A. The list of physical obscrv- 
ables including the above operators both in electronic and 
photonic systems are given in Table I for comparison. 

Finally, it should be noted that the optical Hall ef- 
fect comes from the particle-wave duality of an optical 
wavepacket and the geometrical/topological property of 
a wavcfunction. Therefore, this effect can be observed 
in a macroscopic wavepacket of light described by classi- 
cal electrodynamics, when a wavepacket under consider- 
ation is approximately coherent. In this sense, the sec- 
ond quantization is not always necessary. The second 
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TABLE I: Operators relevant to wavepacket dynamics. 



Electronic system 



Photonic system 



J dr rip^ {r)tl){r) 
J dr ip''{r) [—iVr — eA{r) + eB x r] ■0(r) 
/ dr il)^(r) [r x (— iVr) + s] tl){r) 
J dr {r)ev{r)ip{r) 
f dr tp''{r) [f r X v{r) + g/iss] i^ir) 
i / rfr '(/'t(^)|ij(j.),r|v(r) 
f /dr [ff,V'"''(r) |i?(r),r} V^(r)] 
^Jdrrx [i /,V>t(r)|H(r),r|^(r)] 



i / dr [E{r) ■ D{r) + H{r) ■ B{r)] 

undefined 

i / dr [D{r) x B(r) - B(r) x D{r)] 
^Jdrrx [D{r) x B(r) - S(r) x D(r)l 
undefined 

undefined 

^Jdrr [E{r) ■ D{r) + H{r) ■ B{r)] 
i / dr [E{r) X if (r) - H{r) x i;(r)] 
^Jdrrx [E{r) x (r) - (r) x E{r)] 



H 
R 
P 
J 
I 

M 

n 

V 

J 



H: Hamiltonian, H{r): first-quantized Hamiltonian 

R: position (dipolo moment), P: (pscudo) momentum, J: angular momentum, s: spin matrix 

I: charge current. M: magnetic moment, v{r): first-quantized velocity 

TZ: position weighted l)y ('ii(>] g>' density. 7-": energy cuirt'iit. JT: rotation of energy current 



quantization is adopted, for convenience, to calculate a 
motion of a wavepacket on an equal footing with that of 
an electronic system, as shown in Sec. II C. As long as we 
consider an approximately coherent wavepacket in a sin- 
gle particle approximation of quantum theory of photon 
or in a linear approximation of classical electrodynamics, 
results obtained by both formalisms coincide with each 
other as shown in Sec. Ill A. Detailed remarks on the 
relation between quantum and classical pictures of the 
optical Hall effect is given in Appendix A. 



B. Eigenfunctions in a periodic system 



Here we introduce eigenfunctions in a periodic system, 



ih-r—iE-nut 



(20) 



where F = E or H. The symbols n, A and k represent the 
band index, the index for degenerate modes in the n-th 
band and the lattice momentum, respectively, and Enk 
is the energy eigenvalue of the n-th band, which may be 
degenerate. It should be noted that the band index n is 
not needed in locally isotropic systems without periodic 
structure, but we must keep the index A to distinguish 
different polarization states. U^^^{r) and U^^^{r) are 
Bloch functions for electric field and magnetic field, re- 
spectively. It should be noted that the lattice momentum 
fc will be restricted to the first Brillouin zone in the rest 
of this paper. The eigenfunctions satisfy the Maxwell 
equations. 



M(r)at*5,fc(r,f) 



V.r X * 



H 

iXk 



(21a) 
(21b) 



and they lead to the following eigen equations, 

V. X [Jl-'V X ^^^kir)] = e(r)£^fe*f,fe(r)(22a) 

V. X [e-^V X ^^M] = nr)El„^^Mm 

In the case of *e^{r) = *e{r) and *jj^{r) = /u(r), we can 
orthonormalize the Bloch functions with the same lattice 

momentum k as. 



i 



dr 



ws ^^ws 
dr 



ws '-"ws 



f/nAfe(r-)e(r-)C/„^vfe(r-) = <5„„'(5AA',(23a) 
U^xlir)riir)U^,ykir) = 5w^AA',(23b) 



where the domain of integration is the unit cell with the 
volume vws ■ The orthonormality for the eigen functions 
will be discussed later. 

We introduce the Fourier transformation, 



U^XkiG) = / 
jy 



dr 



-iG-TTTF 



ws ^nvs 
dr 



l/„^fe(r), (24a) 



i(G-G')-r 



WS "WS 



e(r), (24b) 



m(G,G') = / ^e-(«-«')"-M(r), (24c) 
Jws 

where F = E or H. G represents a reciprocal lattice vec- 
tor. In terms of the above representations in the Fourier 
space, we introduce the following compact representation 
for the latter convenience, 

\U) = [U{Go),UiGi),U{G2),---], (25) 
and the tensors, 



Pk{G, G') = {k + G)5{G, G') = K5{G, G'), 



[Si]^, {G,G')^-ie,,kS(G,G'), 

sf = • sri-'Pk ■ s, 

=Pk- S7-'Pk ■ S, 



(26a) 
(26b) 

(26c) 
(26d) 
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where we have introduced the abbreviation K ^ k -\- G. 
The inner product of the Bloch functions and the algebra 
of the above tensors are represented as, 

dr 



where the fe-integration is over the first Brillouin zone, 
i.e., 



{U\V) 



{U\iS\V) 
■iPk-S\U) 



J-9 



ws '-"ws 
dr 



-U*{r)-Vir), 
U*{r) X V{r), 



/ws 

[KxU{Go),KxU{Gi), 
Thus the orthonormaUty is rewritten as. 



(27a) 

(27b) 
■](27c) 

(28a) 
(28b) 



dk = [ 

BZ Jfcelst BZ (.^""j 



(35) 



The operators a„Afe and al^^f^ are defined by 



dnXk 



By this notation convention, the Maxwell equations 
for the Bloch functions are represented in the following 
compact forms. 



(KmU^Xk) 



nXhi 1 



iPk 



-iPk ■ S\U, 
{Km^xk)=0, 



nXkl ' 



where 



\K) = [O,--- ,0,fe + G,0,---]. 



(29a) 
(29b) 
(29c) 

(30) 



P^(CA/e|e|£;) = T^(*fAfcl/^l^f) 

I dr [*f,\(r) ■ D{r) + *^4(r) ■ B{r)] , 

(36a) 

^nk ^nk 

= jdr [D{r) ■ ^^^^{r) + B{r) ■ *f,;,(r)] , 

(36b) 

By using Eq. (17), the following commutation relation is 
obtained. 



^nXk 



O-nXk, C^n'X'k 



6nn'SxX'S{k - k'). (37) 



From Eqs. (29a) and (29b), we can easily derive the fol- 
lowing equations, 

^kP^Xk) = ^^KklU^Xk), (31a) 

^k\U^xk) = T^ElklU^xk)- (31b) 

In rclativistic systems, the orthonormality for the 
eigenfunctions are conventionally represented in terms of 
the inner product defined by 

{f\g) = I J dr[r{r,t)-{dtg{r,t)]-[dtrir,t)]-g{r,t)], 

(32) 

and we obtain the following orthonormality relation as 
shown in Appendix B, 

i^^Xk\*^\^n'X'k') = Snn'Sxx'Hk ' k'), (33a) 
i^^Xkm'^n'X'k') = Sun'SxX'Hk ' k'), (33b) 

i^txkm^^'X'k') = (*^4ImI*?a'/c') = 0. (33c) 

where S{k-k') = {2Tr)^S{k-k'). It should be noted that 
we have used *e'^{r) = *e(r) and *jl'^{r) = Jl{r) in the 
derivation of the above relations. This orthonormality 
is required to expand the electric and magnetic fields in 
terms of the eigenfunctions as follows, 



^nxkir,t)al^kh4a) 



C. Equations of motion 

In order to see the effect of geometrical phase on the 
trajectory of a wavepacket, a driving force is needed [53]. 
This is because the geometrical effect is given by the vec- 
tor product between the driving force and the Berry cur- 
vature as we shall see later. In an electronic system, 
a driving force is most conventionally produced by the 
gradient of electric potential. The counterpart in a pho- 
tonic system is given by the gradient of *e(r) or *jl{r). 
Of course, a periodic structure itself gives a gradient. 
However, this effect is exactly taken into account by con- 
sidering optical Bloch states as shown in Appendix C 
where readers can find details about the basic features 
of an optical wavepacket in a periodic system. Thus we 
regard the deviation from a periodic structure as a driv- 
ing force for the Bloch states and treat it perturbatively. 
Here, the perturbation is introduced as a modulation su- 
perimposed onto the periodic structure of *e (r) and (r) 
as 



e-'(r-)-%'(r-)e-i(r), r'(r) 



E{r,t) = V / dkEr^k 
n,x -^Bz 

X [^nXk{r,t)anXk- 

H(r,t) = V / dkEr^k 



(38) 



n,A 

X 



where 7e(r) and 7^(r) are scalar functions, and we call 
them "the modulation functions" hereafter. This kind of 

modulation does not change the local symmetries of e (r) 
and /I(r), and does not violate the energy conservation. 
We summarize the definitions of Berry connection and 
curvature in Table II and the main results obtained here, 
i.e., the equations of motion for an optical wavepacket. 



7 



TABLE II: Berry connection and curvature. 



Electronic system Photonic system 



Bloch function \Ur,xk) Wi^xk), W^^xk) 

Normalization {U^xk\Ur,' x'k) = <5„„"5aa' {U^xh\*^\U!i^,^,^) = {U^xkW^W^^xk) = 3„r,'Sx 

Berry connection [Ank]xx' = -i{U„xk\'V hU„x' k) A„fc = i [Af^. + Af^,] 

' ^' i{U^xkmVkU^x'k) 



^nkl 



XX' 



Berry curvature n„fc = V*. x A„k + jA„fe x A„fc il„h = V*. x A„fe + iA„fe x A„fe 



in Table III. Appendix D supplements details about the 
commutation relations and expectation values of various 
operators which are needed to derive the equations of 
motion in a modulated system. 

Now we derive the equations of motion for the dynam- 
ics of an optical wavepacket. An exact wavef unction \^) 
satisfies the Schrodinger equation, 



urn- 



This equation is derived by applying the variational prin- 
ciple to the quantity [32], 



(40) 



It is natural to consider that the trajectory of the 

wavepacket is determined in terms of the effective La- 
grangian which is given by replacing with a varia- 
tional wavepacket, \W), characterized only by the cen- 
ters Tc, fee of the position and wavevector, respectively, 
and the polarization state [33]. Although \W) can be 
brought closer to j'l') by enlarging the number of vari- 
ational parameters, those concerning the details of the 
wavepacket are neglected here. This approximation is 
justified when a modulation is weak and slowly varying, 
and has been successfully applied to the quantum Hall 
system and gives the semiclassical understanding of the 
motion of magnetic Bloch states [34, 35, 54]. 

In general, the modulation may mix the creation and 
annihilation operators. Here we consider the situation 
in which the modulation is sufficiently weak and time- 
independent, and this mixing is negligible. Therefore, 



J 



the approximated wavepacket can be constructed as 

\W) = / dkw{k,kc,rc,Zc,t)y2^cxalxk-rJ^)' 
Jbz , 



w 



(fe, fee, r„ Zc, t) = Wr{k - feJe-^'^C'''--^-*), (41a) 



where Wr{k—kc) is a real function, and Wr{k—kc) and Zc\ 



(39) satisfy the normalization conditions, /^^ dk (k — kc) 



1 and I^caI^ ~ 1, respectively. We assume Wr{k — 
kc) has a sharp peak around fee = dk w^{k — kc)k. 
Here we require that the center of wavepacket, Tc, is self- 
consistently determined by 

Tc = I dkwl{k-kc)\S/k'&{k,rc,Zc,t)-{zc\^^„k\zc)\, 
Jbz 

(42) 

where A„fe is the Berry connection defined by 



A„fe = - [A^fe + A^fe] , 

[aSc]aa' = -i^U^xkmVkU^x,^), 
and we introduced the abbreviation, 

{Zo\M\Zc) = ^Z*aMaa'^;cA'- 
A,A' 



(43a) 

(43b) 
(43c) 



(44) 



The annihilation and creation operators of approximate 
eigen modes are defined by 



flnAfe;ro — 



a, 



n\k\rc 




(45a) 
(45b) 



r 



These operators satisfy the same commutation relation as that in a periodic system, 

anAfe;rc5 <^l'A'fc';r, 



5nn'5xx'hk - fe')- (46) 
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The approximate eigen modes depend on the variable 
Tg. Thus we must consider also the operator ^ro^^nxk-r 
when estimating the effective Lagrangian. However, we 



can show that the contribution from Vi.„a, 



by using 



a. 



rc"'n\k;r. 

= and a„Afc;rjO) 



vanishes, 
0. 



As was mentioned above, we consider the situation in 
which the mixing of the approximate annihilation and 
creation operators is negligible and assume that the re- 
lation anXk;rjO) = 0, holds. The expectation values of 
H and TZ are estimated by the derivative expansion with 
respect to 'ye{r) and 7/i(r) as 

{W\H\W) S £r,k^,r,,,^ (47a) 

{W\n\W) ^ Enk,;rA^kJ{kc,rc,Zc,t) - {Zc\A„kJZc)] , 

(47b) 



where Enk;rc = le{rc)lu.{rc)Enk, 



£nkc;rc',z 
EnkctTc 



Vr, In 



(^c|A„feJ^c),(48) 



and 



Anfe = ^ [Affe - A, 



H 



This function A„fc is a difference between the Berry con- 
nections of the electric and magnetic parts. Although 



(^c|A 



E,H\ 



,) depends on the representations of eigen 



modes, {zc\^nk\zc) does not. This issue is related to the 
gauge transformation in fc-space given in Appendix C2. 
The property of {zc\^nk\zc) is also discussed in Ap- 
pendix E. 

In the above evaluations, we assumed that the shape 
of w'^{k — kc) is sufficiently sharp compared to the slow 
variations of Enk-^r^ and A„fe(fc) around fc^, and neglected 
terms which depend on the shape of u'^,(k — fee). There- 
fore, even with the pcrturbative modulation, we can re- 
gard Tc defined by Eq. (42) as the center of gravity 
(M^|7?.|M^)/(VF|iI|H^). In the present approximation, we 
assume that the modulation is so weak and smooth that 
we can neglect the second order derivatives of the modu- 
lation functions in the effective Lagrangian and the equa- 
tions of motion, which we shall derive later. In the equa- 
tions of motion, the difference between r,, and the center 
of gravity appears as higher derivatives than the original 
derivatives. Therefore we may neglect the difference be- 
tween Tc and the center of gravity due to the derivatives 
of the modulation functions, while we cannot neglect the 
correction due to the first derivatives in 

Here we introduce the effective Lagrangian in order to 
derive the equations of motion of a wavepacket. 



i^eff = {W\ij^-H\W). 



(50) 



The first term on the right-hand side of Eq. (50) is cal- 



culated as follows. 



iWVj,\W) 



kc -Tc- kc - {Zc\AnkJZc) + i{Zc\Zc) 

d 



+ 



dt 



Jbz 



(51) 



Neglecting the total time-derivative, we obtain the final 
form of the effective Lagrangian, 

-^efF — * * (•^■cl Ayife^jZc) "l" i(^Zc\Zc) £nkc;rc;Zc- 

(52) 

From this Lagrangian, the equations of motion are de- 
rived as follows, 



'f'c — ^ kc^nkc;rc:Za + fee X (^c| ^^nfcc I -^^c) 

-«(2c| [f^ ■ A„fe^, A„fcJ \zc), (53 



a) 
(53b) 



14) 



kc ■ Ank, + fc ■ A„fe, 



\zc), (53c) 



(49) 

where 



Vr. In 



leiVc) 



and Ctnk is the Berry curvature defined by 

^nk = Vfc X A„fe + iAnk X A„fc. 



(54) 



(55) 



It should be noted that the above equations of motion 
satisfy the energy conservation, i.e., 



dt 



(56) 



Rigorously speaking, a Lagrange multiplier is needed for 
the constraint {zc\zc) = J2x l-^cAp = 1 in the derivation 
of the equations of motion, while the constraint is implic- 
itly imposed here. Therefore, in the above equations of 
motion, we should consider that the above constraint is 
always imposed. 

In generic cases with periodic structures, we cannot 
analytically evaluate A„fc, ^Ink nor A„fc. However, in 
principle, we can numerically calculate them and thus 
solve the equations of motion of a wavepacket subject 
to a modulation superimposed onto a periodic structure. 
Appendix E presents formulae for ilnk and A„fc as well 
as the internal rotation of a wavepacket, which are use- 
ful for numerical calculations. Although the relation, 
{zc\A^^\zc) = (zc|A^^|zc), is not generally proved, this 
can be confirmed at least for some systems with locally 
isotropic e(r) and iJ,{r), e.g., for the elliptically polarized 
light in systems without periodic structure and for the 
TM modes in two-dimensional photonic crystals where 
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/X is constant. In this case with A„fc = 0, we can re- 
place the perturbed energy as £nkc;ra\zo ~^ Enk^^r^ — 

Finally, it is noted that the optical Hall effect is orig- 
inated by the second term on the right-hand side of 
Eq. (53a), which is sometimes called "anomalous veloc- 
ity" . The anomalous velocity is the vector product of the 
Berry curvature (zdJ^nfe^ |zc) and the driving force fee, 
i.e., the gradient of a superimposed modulation. There- 
fore, both of the driving force and the Berry curvature 
are needed for this phenomena. In an electronic system 
under a strong magnetic field, it is pointed out that the 
Berry curvature is closely related to the internal rotation 
of an electronic wavepacket [34]. Thus, a similar relation 
is also expected in a photonic system. Actually, in Ap- 
pendices C and E, we can see the close relation between 
the Berry curvature and internal rotation of an optical 
wavepacket. 



III. APPLICATIONS 

A. Transverse shift in reflection and refraction 

As an application of the theoretical framework devel- 
oped in Sec. II, we consider the case with locally isotropic 
e(r) and ^{r). In this case, we can write down the equa- 
tions of motion (53a)-(53c) in simple forms. For this 
purpose, we first calculate the gauge field Afe and the 
Berry curvature ilfc. In a helicity basis, the eigenvectors 
for the right(+)- and left(— )-circular polarizations can 
be written as 



1 



(egiie^), U. 



H 
±.k 



1 



effect induced by a modulation of refractive index [20]. 
The equation for \zc), Eq. (60c), describes a phase shift 
by the directional change of propagation discussed in 
Refs. [17-19]. This equation gives the solution |z°"') = 



2+ , 



], where [zj,") — 



+ ' ■ 



(e^^ieg), (57) 



of polarization. Q is a solid angle made by the trajec- 
tory of momentum: Q — § dk- [Afc]+-|_ — dSk ■ [$7^]++ 
where dSk is the surface element in fc-space and 5* is a 
surface surrounded by the trajectory. Our approach can 
be easily generalized to treat systems with periodic struc- 
tures on the same footing, and offers a powerful tool for 
applications compared with the eikonal approximation 
[30]. 

The simplest example of the optical Hall effect is re- 
alized as the transverse shift at the interface refraction 
and reflection. There have been a number of studies on 
the shifts within and out of the incident plane at the 
total reflection. The former is well known as the Goos- 
Hanchen effect [55, 56] and has been explained in terms 
of the evanescent wave penetrating into the forbidden re- 
gion. The latter one, which is referred to as the Imbert- 
Fedorov shift, was interpreted by F. I. Fedorov [21] as an 
analog of the Goos-Hanchen effect and was observed ex- 
perimentally by C. Imbert using multiple total reflections 
[22] , followed by a number of theoretical approaches [23- 
25]. Furthermore, it was pointed out that the shift out 
of the incident plane could also occur in partial reflection 
and refraction [25-27]. However, some of the theoretical 
predictions for the amount of shift contradict each other. 
One reason is an experimental difficulty for a measure- 
ment of the tiny shift, as the shift is only a fraction of a 
wavelength. It was only recent that the Imbert-Fedorov 
shift is measured for a single total reflection [28]. Thus 
the physical mechanism for the transverse shift is still 
controversial. 



where eg^^ with — k/k are the orthogonal unit vectors 
in the spherical coordinate of the fc-space. After some 
calculations we obtain 



A^ - A« - 



(58) 



which then yields [57] 
cos 9 



fcsin^"^^^' = (59) 



Thus the equations of motion are simplified as follows, 
fc 

r-c = v{rc)-r^ + kc X {zc\flkjzc), (60a) 



kc = [V^^v(rc)] kc, 
\zc) = -ike ■ Akjzc), 



(60b) 
(60c) 



where v{r) — 1 / ^/eirjjjjr) , and \z) — [z+,z_] is rep- 
resented in a helicity basis. As is analogous to the 
Hall effect in electronic systems, the second term on the 
right-hand side of Eq. (60a) describes the optical Hall 




FIG. 1: Transverse shift of light beams in the refraction and 
reflection at an interface. 

In our previous paper [20], we calculated this trans- 
verse shift by using the conservation of the z-component 
of TAM for individual photons, which follows from the 
equations of motion (53a)- (53c) applied to this interface 
problem. Bliokh et al. [50] then questioned the result, 
claiming that it does not match their result of the shift 
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TABLE III: Equations of motion of optical wavepacket. 

Modulation e"^(r) -» 7?(r)e~^(r), /i"^(r) -» 72(r)/u"^(r) 

Effective Lagrangian Lea = kc ■ r-c — kc ■ {zc\A„k^\zc) + i{zc\zc) — £nk^;r^;z^ 

Equations of motion rc = Vfc„fnfcc;T-c;zc + x {zc\ilnka\zc) — iizc\ [ft ■ A„fc^,A„fc^] \zc) 

kc — ^ rc^nkclVc'jZc 

\ic) = -like - A„fe^ + ft ■ A„fej2c) 



Perturbed energy 






[-^nka 


' ^nkr} 


2c)J Enkcirc Ejikc'.rc ^^ fc 


(Zc\A„kc\Zc) 




Enkc'-Tc — 


= je{rc)'yij.{rc)Enkc, ft = - 


[Vrc In 


Itir-c) 1 


1 EnkciTc, ^nkc = 1 [-^nko ~ 


Anfec] 


Energy conservation 


d c 


=0 











for elliptically polarized Gaussian beams. Here we show 
that, in all cases with generic polarizations, an identical 
result for the transverse shift of each beam is given by the 
following different approaches, (i) analytic evaluation of 
wavepacket dynamics, (ii) TAM conservation for individ- 
ual photons in Ref . [20] , and (iii) numerically exact simu- 
lation of wavepacket dynamics. It agrees with a result by 
classical electrodynamics, as presented in Appendix F. In 
Sec. IIIB, we shall resolve the inconsistency between the 
identical result by the approaches (i)-(iii) and that given 
in Ref. [50] . Thereby the validity of our theory presented 
here is completely guaranteed. 

Our equations of motion are not directly applicable 
to the refraction and reflection problem at a sharp in- 
terface since they require the slowly varying conditions 
[Vlnej, |Vln/i| <C k. Indeed, our equations of mo- 
tion do not correctly describe a splitting of an incident 
wavepacket into reflected and transmitted wavepackets 
nor changes of their polarization states at the interface. 
However, in a case with a flat interface, a simple exten- 
sion of our theory works well, as will be explained in 
the following. Here, as shown in Fig. 1, we consider the 
case where the incident beam comes from the region with 
X < and 2; < along the plane of y = const., and the 
interface is the z = plane. 



1. analytic evaluation 
Far from the interface, Tc in Eq. (42) is easily estimated 

as 

where A = I for t —00, A = T ov R for t 00. The 
momenta k^''^'^ and the polarization states \z^''^'^) are 
those of the incident (/), transmitted (T) and reflected 
(i?) beams, respectively. Due to the wavepacket splitting 
at the interface, our semiclassical equations of motion do 
not tell us the values of \z"^) and 1}^. Hence, for these 
variables, we borrow the results of reflection/refraction 
of a plane wave. With a natural choice of a wavepacket 
presented in Appendix F, the t/-componcnt of the first 
term in Eq. (61) is unchanged at the interface. Thus, the 



transverse shift comes only from the second term as 

6y^ = -(2^|Afe.|z^) + (z^|Afe.|/). (62) 

where A = T or R. Substituting Eq. (59) in each Berry 
connection {z^\AkA\z^) {A = I,T,R), we obtain the 
following equation for the transverse shift 

= ^7^[(^'*|a3|^^)cos0A-(z^|a3|zOcos^?,], 

(63) 

where A = T 01 R, Oj t.b. arc the angles between the 
positive z-axis and the propagating directions of the in- 
cident, transmitted and reflected beams, respectively. 

2. total angular momentum conservation 

Then, what is the physical meaning of the above result? 
Firstly, it should be noted that (z|(T3|z) = |z+|^ — [^-l^ 
represents the magnitude of spin polarization in the di- 
rection of fee, i.e., (z|(T3|z) = ±1 for right/left-circular 
polarizations, (z|cr3|z) = for linear polarizations, and 
|(2:|a3|2;)| < 1 for elliptic polarizations. Therefore, it is 
intuitively expected that this phenomena is closely re- 
lated to the angular momentum of wavepacket. For a 
system with rotational symmetry around the z-axis, the 
equations of motion lead to the conservation of the z- 
component of the following TAM, 

k 

jc = rc X kc + {zc\a3\zc)-rr- (64) 

This conservation is expected to hold even in the case of 
a sharp interface, because it is based on the rotational 
symmetry around the z-axis. Actually, from Eq. (63) for 
the transverse shift and Eq. (64) for the TAM, we can 
reach the conservation of the z-component of the TAM 
for each of individual photons 

ji=f.,ji=3z, (65) 

where j^'^'^ are the TAM of incident, transmitted and 

reflected beams, respectively. It just makes sense that 
the incident beam is regarded as a collection of photons; 
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each photon is reflected or transmitted stochastically at 
the interface. As we shall see in Sec. IIIB (and Appendix 
F in detail), Eq. (63) for the transverse shift is consistent 
with the result derived in classical electrodynamics. In 
this sense, this photon picture is implicitly incorporated 
already in classical electrodynamics. 

Inversely, assuming the conservation of TAM for in- 
dividual photons, we can derive the transverse shift as 
Eq. (63). This is what we have done in our previous pa- 
per [20] . This derivation of the transverse shift is akin to 
the derivation of Snell's law based on the particle picture 
of light in which the refracted and reflected angles is ob- 
tained from the conservation of energy and momentum 
(parallel to the interface) for individual photons. 



3. numerical simulation 

In order to verify our theory quantitatively, we check 
the property of the transverse shift in more detail. To 
obtain (^i^lcrsl^;"^), we decompose the incident wave as 



1^ ) = r- \P) + +^ \S), 



V2 



(66) 



where \p) = -^[1, 1] and \s) — -!=[— represent the p- 
and s-polarized states. Straightforward calculation yields 

(2^1^31^^) 

_ 2 [{z'\a3WmA;As) + {z'\a2W)'^{AlA,)\ 

- [1 + {ZI\CT,W)]W + [1-{ZI\CT,W)]W 

(67) 

with ^ = T or i?, and Tp and Tg {Rp and Rs) are the 
amplitude transmission (reflection) coeflacients for p- and 
s-polarization, respectively. 

When we focus on the partial reflection and refraction, 
Ap and Ag are real, and Eq. (63) is rewritten as follows. 



tan 9 1 

2ApAs cos 9a I cos 9i 
{\ + {z^\o^z^)\Al^{\-{z^\a^z^)\Al 



- 1 



qs) 



where A = R (yvT . This means that the incident beams 
with \z) = [e'F*^ cos |, e^*2 sin |], where 9 and </> rep- 
resent the spherical coordinate of the Poincare sphere, 
cause the shift of the same magnitude and the same di- 
rection with each other, i.e., the shift independent of 
the sign of In addition, the incident beams with 
— sin |, e^'^ cos |] cause the shift of the same 
magnitude as the above beams, but of the opposite di- 
rection to them. In the partial reflection and refraction, 
no shift is observed for the incident beam with linear 
polarization. 

On the other hand, in the total reflection, we have 



represented by 

= - 



1 



tan 9 1 



[(zV3k')[3?(i?;i?s) + i] 

+ (zV2|z^)3(i?;i?.)]. (69) 



\Rv 



\Rs 



1, and the shift of the reflected beam is 



In particular, for the incident beam with linear polariza- 

tion \z) = ) ^^]' shift is the same magnitude 

and the same direction for ^ = a and (j) = n — a. The 
direction is reversed by the replacement (p —<j) without 
change of the magnitude. 

We have confirmed all the above features quanti- 
tatively by numerically solving Maxwell equations for 
wavepackets. In Ref. [20], we have presented the results 
only for the incident beam with right- circular polariza- 
tion. Here, to complete the argument, we present the re- 
sults of the numerical simulations for more generic cases. 
Figure 2 shows the shifts for the incident beams with the 
elliptical polarization z^j^/zL = 2 at the interfaces with 
relative refractive indices (a) n = 2.0, (b) n = 0.8, and (c) 
for the incident beam with linear polarization z^j^/ z[_ = i 
with n = 0.5. (We take the value of magnetic permeabil- 
ity common in both media upper and lower the interface, 
i.e., fii = /i2, in these simulations.) The solid and dashed 
lines represent the analytic results Eq. (63) for transmit- 
ted and reflected beams, respectively. The filled circles 
and squares are the results of simulations for transmit- 
ted and reflected beams. We note that, in Fig. 2(c), the 
shift for the linearly polarized beam is nonzero only for 
a region of total reflection, in accordance with our ana- 
lytic result. In all cases, the numerical results excellently 
agree with Eq. (63), thus verifying our theory. (We have 
confirmed this consistency also in cases where both of 
permittivity and permeability are different in two media 
upper and lower the interface, i.e., ei ^ £2 and fii ^ /i2.) 

Finally we should comment on our constitution 
method of a set of incident, transmitted and reflected 
wavepackets, which is an exact solution to Maxwell 
equations. In each numerical simulation, we have con- 
structed an elliptically-polarized incident wavepacket as 
a superposition of plane waves with a common polar- 
ization state, i.e., \z^) which is independent of k. This 
is a natural definition of incident wavepacket. Oth- 
erwise, the concept "an elliptically-polarized incident 
wavepacket" gets fuzzy, and the linear composition from 
and decomposition to different orthonormal bases of in- 
cident wavepackets are violated. This is because inci- 
dent wavepackets with different functions of {|2;^(fe))}'s 
for constituent plane waves can have the same mean po- 
larization state \z^{k^)). Imposing the exact boundary 
conditions, transmitted and reflected wavepackets are au- 
tomatically generated. For partial reflection, a single 
incident wavepacket split into reflected and transmitted 
wavepackets after reflection/refraction at the interface. 
The position of each wavepacket is estimated when each 
wavepacket is far from the interface. It should be noted 
that the numerical simulations exactly take into account 
the changes of shapes of wavepackets, while the analytic 
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FIG. 2: Shifts of reflected and transmitted beams, n is the relative refractive index of the upper medium with respect to the 
lower medium. Xi is the wave length of incident light in the lower medium. The solid and dashed lines represent the analytic 
results Eq. (63) for transmitted and reflected beams respectively. The filled circles and squares are the results of simulations 
for transmitted and reflected beams. 



evaluation assumes the sharpness of a weight function for 
the superposition. 



dition for a set of wavepackets in paraxial approximation. 
This boundary condition is different from the correct one: 



B. Remarks on other theories 

Recently, Bliokh et al. calculated the shift for an elhp- 
tic Gaussian incident beam in classical electrodynamics, 
and their result disagrees with that obtained from our 
theory [50] . They attributed the difference to a "fallacy" 
in our TAM conservation for individual photons. We ex- 
plain below in detail that our theory is totally free from 
the criticism. To prove this, it is enough to show that 
Eq. (63) is equivalent to the transverse shift evaluated 
in classical electrodynamics, i.e., the result by Fedoseev 
[26, 27]. His procedure of calculation is as follows. First 
construct a wavepacket by a linear superposition of plane 
waves. By taking into account the exact boundary con- 
ditions of electromagnetic fields (Eqs. (70a) and (70b)) 
at a flat interface, as in the textbooks of optics or clas- 
sical electrodynamics [38, 56], one can construct trans- 
mitted and reflected wavepackets as an exact solution 
to Maxwell equations. The center of each wavepacket is 
defined as an average position weighted by each energy 
density. The result of this calculation by classical elec- 
trodynamics [26, 27] is identical with that by our the- 
ory (Eqs. (68) and (69) derived from Eq. (63)) in the 
second quantized formalism. The details are presented 
in Appendix F. Hence we checked that the following 
three approaches give the same transverse shift for each 
wavepacket; (i) analytic evaluation of wavepacket dynam- 
ics both in classical and quantum-mechanical formalisms, 

(ii) TAM conservation for individual photons (Eq. (65)), 

(iii) numerically exact simulation of wavepacket dynam- 
ics. 

There remains an inconsistency between the identical 

result obtained by (i)-(iii) and one by Bliokh et al. in 
Ref. [50]. One reason is an inappropriate boundary con- 



t ■ [E\r, t) + E^ir, t)\=t- E'^{r, t), (70a) 
t ■ [H\r, t) + H\r, t)\=t- H^{r, t), (70b) 

where t is an arbitrary unit vector parallel to the inter- 
face, E"^ and {A = I,T,R) are electric and magnetic 
fields of incident (/), transmitted (T) and reflected (i?) 
beams, respectively. Another reason for this contradic- 
tion comes from the deflnition of a center of wavepacket 
in Ref. [50]. The methods (i)-(iii) commonly use the 
position averaged with a weight of energy density. In 
Ref. [50], on the other hand, the center is defined as a 
center of the wavepacket projected onto its mean polar- 
ization state. The center of the wavepacket in the for- 
mer definition, i.e. the position averaged by the energy 
density, can be easily measured by photon counting, as 
employed in two measurements on the Imbcrt-Fedorov 
shift [22, 28], while the latter deflnition requires count- 
ing of photons projected onto a specifled polarization. 
The agreement between totally different approaches, i.e., 
(i)-(iii) and classical electrodynamics, suggests that our 
deflnition is a natural one. 

Finally we should comment on the relation between the 

conservation laws of TAM in the wave and particle pic- 
tures of light. In Ref. [50], it is claimed that, for an inci- 
dent beam with an elliptic; polarization, the conservation 
of TAM for individual photons (Eq. (65)) is inconsistent 
with the conservation of TAM for whole beams, 



•^^ niiJL2 cos e^-'"' 



(71) 



However, as we shall show below, Eq. (65) is a sufficient 
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condition for Eq. (71). From Fresnel formulae, we have 

^2 n2^J.l cos 9t 



^ ' ni/i2 cos6'/ ' 

-, r,2 , ^2 COS0T 

i = Kg + Ig , 

ni^2 cos di 

for the p- and s-polarized beams, and also 

^2 ?^2Aii cos 0T 



1 = 



nifj,2 cos 9 1 

for a beam with an arbitrary polarization, where 



l + |mP ' 



l + |mP ' 



(72a) 
(72b) 

(73) 
(74) 



and m = z^/zp, \z) — Zp\p) + Zg\s). The above for- 
mula represents the conservation of energy flow or equiv- 
alently the conservation of the number of photons. One 
can easily see that the above formula and Eq. (65) yields 
Eq. (71). To summarize, the TAM conservation for indi- 
vidual photons (Eq. (65)) has neither contradiction nor 
inconsistency with other theories. 



C. Two-dimensional photonic crystal 

We consider two-dimensional photonic crystals, where 
e*(r) and Jl{r) are scalar variables e(r) and fi{r) which 
are periodically modulated in the xy-plane. Along the 
^-direction e(r) and ^J-{r) are assumed to be uniform. 
However, it is noted that in general there may appear 
ordinary degeneracies at symmetric points and accidental 
degeneracies at some specific points in the Brillouin zone. 
Around these points, the semiclassical argument based 
on the adiabaticity would not be a good approximation, 
and it is needed to seriously incorporate the dynamics of 
the wavepacket. We restrict ourselves to bands without 
degeneracy here for simplicity. In this case, the inversion 
symmetry of the periodic structure must be broken in 
order for a band to have nonzero Berry curvature. This 
is because the Fourier transformation of e(r) and /i(r) are 
real-valued, when a system has the inversion symmetry. 

For a wavepacket constructed from a non-degenerate 
band, the equations of motion in Eqs. (53a)-(53c) are 
reduced to the following ones, 



'f'c — ^ ka^nka;ra + fee ^ ^nkc 
fee ^ rc^nhc'^Vc^ 

kr ■ A„fc„ -I- ■ A„ 



(75a) 
(75b) 

(75c) 



In the above equations of motion, the most important 
and controllable quantity is the second term on the 
right-hand side of Eq. (75a), i.e., the Berry curvature 
^^nfcc- The anomalous velocity fee x r2„fc^ of the optical 
wavepacket leads to the optical Hall effect. Compared 
with this term, the other correction due to A„fc are 



small as shown below (see Appendix H also). The pa- 
rameter Zc becomes a simple complex number and just 
represents a phase shift. Thus an optimal design for the 
enhancement of the optical Hall effect is equivalent to 
the enhancement of the magnitude of the Berry curva- 
ture. In the present case, the Berry curvature comes 
from an interband effect due to a periodic structure with- 
out inversion symmetry, and roughly scales as the inverse 
square of a band splitting (see Eqs. (E3a) and (E3b)). 
Therefore, we can expect the enhancement of the optical 
Hall effect for wavepackets constructed from Bloch waves 
around nearly degenerate points in the Brillouin zone. 
In two-dimensional photonic crystals, Bloch waves prop- 
agating along the xy-plane [kz — 0) are classified into 
the transverse magnetic (TM) and the transverse elec- 
tric (TE) modes. In other words, the Maxwell equations 
for the Bloch functions (29a)- (29c) decouple into two sets 
of equations, one for the TM and the other for the TE 
modes, and the problem of wavepacket dynamics can be 
more simplified. Appendix G gives useful formulae for 
the Berry curvature r2„fc and A„fc in such modes. As 
for other modes and more generic case with degenerate 
bands, we must use formulae given in Appendix E. 



(a) £ix,y) 



(b) Photonic band = 0.5) 




FIG. 3: (a) Dielectric function and 
a two-dimensional photonic crystal, 
shown in Figs. 4 and 5(b) 



(b) band structure of 
The Brillouin zone is 



We present examples of the Berry curvatures and the 
internal rotations of non-degenerate bands in the two- 
dimensional photonic crystal with ^ — and 



-\r) 



3(5 + 12|^H-8C2) 



3 

E 

i=l 



27r, 

C - cos(6j + ^) 



n 2 



,u 27r; 
C + cos(6j "^^ y) 



(76) 



where bi = (2|^, -f ), 62 = (0, f ), = -61-62, and 
a is the lattice constant. It is noted that, for < |^| < 1, 
^ represents the degree of inversion-symmetry breaking. 
The spatial distribution of e(r) and the band structure 
of TM and TE modes are shown in Fig. 3(a) and (b), 
respectively. 
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Figure 4 shows the Berry curvatures and the internal 
rotations of the first and second bands of TM and TE 
modes. The internal rotation of an optical wavepacket is 
definedby (zc|5„fcjzc) = {W\J\W)-rcXEnk,y k^Enk,, 
where (W^I^TIM^) is the total rotation of energy current 
and the second term represents the rotation of the center 
of gravity (see Appendices C, E and G). We can clearly 
see the correlation between them in each band except 
for their relative sign. The relative sign are roughly de- 



termined by a factor SE = (£'tm(te) nk ~ E' 



'tM(TE) mk) 

at nearly degenerate points fc, where n and m represent 
band indices of nearly degenerate bands. This is because 
the Berry curvature and the internal rotation are propor- 
tional to 1/(5 -E^ and 1/SE, respectively (see Appendix G). 
It is expected that this internal rotation is closely related 
to a physical angular momentum. Therefore, these re- 
sults suggest that we can generate photonic modes with 
angular momentum by using photonic crystals without 
inversion symmetry. 

Before considering the motion of wavepackets in this 
photonic crystal, we should comment on Atm nk and 
^TE nk, which give corrections to energy dispersions 
and group velocities of TM and TE modes. Because 
/i(r) = /io in the present case, it follows from Eq. (G8a) 
in Appendix G that AxMnfe = 0. Thus, when a mod- 
ulation is applied only to the dielectric permittivity 
as l/e{r) 7^(r)/e(r), the energy of the TM mode 
is just rescaled by the factor 7e(rc), i.e., -Etm nfec ^ 
^^TM nfc,;re = 7e(''c)i^TM nfc, ■ On the Other hand. Ate nk 
is nonzero. From Eq. (48), additional corrections ap- 
pear in the energy dispersions of TE modes. However, 
as shown in Appendix H, we can see AxEnfc ^ 0.1a. 
Thus these corrections are estimated to be at most a few 
percent as long as the modulation is sufhcicntly weak, 
i.e., |aVr^ ln7(:(rc)| ^ 1. In the similar argument, we 
can also neglect corrections to the group velocities of TE 
modes compared to their anomalous velocities, at least 
in the present photonic crystal. All the details of this 
issue is given in Appendix H. 

Now we consider the motions of wavepackets con- 
structed from TM and TE modes. It is noted that these 
wavepackets are extended in the z-direction, because the 
z-components of their momentum are fixed as fcj = 0. 
From Fig. 4, we can see that Hfc is strongly enhanced 
near the corners of the Brillouin zone. This enhancement 
is interpreted as a two-dimensional cut of the monopole 
structure in an extended space including parameters [9], 
e.g., ^ in the present case. Therefore, we set the initial fee 
near the corners of the Brillouin zone in order to make the 
effect of anomalous velocity prominent. We superimpose 
the following modulation onto the periodic structure as 
l/e{r)^jUr)/<r), 



1 



1 r 
2 



(n + 1) + (n - 1) tanh ■ 



(77) 



where n > represents a relative refractive index multi- 
plied to the periodic structure in the region a; ^ oo, and 
w is the mean width of the modulation. Here we take 




(e) TM 1st 





(g) 




(h) TE 2iitl 




FIG. 4: Berry curvatures (a,b,c,d) and the internal rotations 
(e,f,g,h) of the first and second bands of TM and TE modes 
in the two-dimensional photonic crystal (^ = 0.5). 



n ~ 1.2 and w — 5a which satisfy the condition of weak 
and slowly-varying modulation. The obtained trajecto- 
ries are shown in Fig. 5. It is found that the shift of Tc 
reaches to dozens of times the lattice constant especially 
for the wavepacket constructed from the TE second band. 

Finally we note that, also in more generic systems than 
discussed above, this effect can be enhanced considerably 
by designing crystal structures. The Berry curvature 
around a nearly degenerate point is determined mostly 
by the splitting, 2|mg|, between neighboring bands. Note 
that the sign of nig depends on details of wavefunctions 
around the nearly degenerate point, while its magnitude 
is determined only by the splitting. Suppose that at 
fc = fco another band comes very close in energy to the 
one considered. The Berry curvature around feg is eval- 
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(a) Trajectory of wavepacket 
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FIG. 5: Trajectories of wavepackets in (a) real and (b) mo- 
mentum spaces. The color of each arrow in (b) corresponds 
to that of each line in (a). The momentum-space trajectories 
in the figure are drawn with appropriate shifts from their ac- 
tual ones which are on the line of ky — or the horizontal 
Brillouin-zone boundary. 



uated as 



(v2|fc_feo|2+TO2)3/2' 



(78) 



where u is a nominal velocity around fco- Thus when the 
light traverses near ko, the shift is estimated as 



=sgn[TOgVa:,7(a;c)], (79) 



where k is the minimum value of |fc — feol when k traverses 
near kg. Therefore the shift is larger for smaller \mg\ and 
K. This argument gives an intelligent explanation of a 
relation among relative magnitudes and signs of the sifts 
in Fig. 5. 



IV. DISCUSSION 

We have presented in detail the derivation of the equa- 
tions of motion for an optical wavepacket within the uni- 
tary theory. In our formalism, the equations are derived 
in the same fashion as those of electronic systems, thereby 
the similarities between them are evident. This suggests 
a broad concept, i.e., topological Hall effect driven by 
the geometrical mechanism, which ranges over a wide 
area of physics such as electronic, acoustic, hydrodynam- 
ical, relativistic and photonic phenomena. For example, 
this optical Hall effect is referred to as "the optical Mag- 
nus effect" [30, 31], in analogy with the Magnus effect, 
which is a transverse aerodynamic effect on rotating ob- 
jects. On the other hand, the relation of this effect to the 
geometrical effect on a spinning particle in general rela- 
tivity was recently pointed out [58]. They derived the 
similar equations of motion to ours by considering the 
motion of a spinning particle in the space with a metric 



gijir) = n{r)5ij. Their argument reminds us that, in 
the early stage of the study on general relativity, Ein- 
stein had tried to formulate the theory by generalizing 
the speed of hght in vacuum. This issue might be related 
to the deep question of the dual nature between the force 
and the velocity in the dynamics. 

Here we should mention the effects called the photonic 
Hall effect [39, 40] and the magnetically induced deflec- 
tion due to the Pitaefskii magnetization [41-46], both of 
which are observed in Faraday-active media subject to 
external magnetic fields. The former effect takes place in 
a random medium, and is theoretically interpreted by the 
magnetically induced off-diagonal components of a diffu- 
sion tensor [39] and experimentally proved to be due to 
the magnetically induced changes in the optical proper- 
ties of scatterers [40]. The latter effect is observed in 
a homogeneous medium [42], and is interpreted by the 
magnetically induced change in the dispersion relation of 
each mode due to the Pitaevskii magnetization [41]. (Ad- 
ditional remarks on these effects are given in Appendix I.) 
On the other hand, the optical Hall effect is caused by 
the anomalous velocity due to the geometrical propriety 
of a wavepacket, which appears without external mag- 
netic field nor scatterers. 

As focused in Appendix E, there is a close relation be- 
tween the Berry curvature and the internal rotation. It is 
physically expected that an internal rotation can be re- 
lated to an internal angular momentum. From this view- 
point, the various modes of Laguerre-Gauss beam, which 
have internal orbital angular momenta [59] , are of a par- 
ticular interest. The Imbert-Fedorov effect is expected 
for these modes as well as for circularly polarized states, 
and theoretical and experimental investigations on this 
problem has been done recently [60-62]. In Sec. IIIC, 
we have shown that there appear photonic modes with 
internal rotations in a two-dimensional photonic crystal 
without inversion symmetry. Angular momentum corre- 
sponding to this kind of internal rotation would be de- 
tected by measuring a torque taken by a photonic crystal 
when we inject a linearly polarized light into the crystal 
(through a buffer layer if needed). In addition, when 
a photonic crystal is composed of Faraday-active media 
and subject to an external magnetic field, there would 
take place the magnetically induced deflection due to the 
Pitaevskii magnetization caused by this kind of generic 
internal rotation, rather than by the spin of circular po- 
larization. 

Lastly, we make a remark on the relevance of the 
quantum nature in the geometrical/topological proper- 
ties discussed in this paper. Although we have formu- 
lated the theory of an optical wavepacket in the quantum- 
mechanical formalism in order to clarify its connection to 
that of an electronic wavepacket, the phenomena itself is 
based on duality between real-space coordinates and mo- 
menta which is common in wave dynamics. Therefore, 
the topological Hall effect is generic in both quantum me- 
chanics and classical wave dynamics. (See the argument 
in Appendix A.) Actually, we can extend the argument 
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presented here to other kinds of wavepackct dynamics, 
e.g., dynamics of sonic wavopacket which is mentioned 
in Ref. [30]. The sonic HaU effect in phononic crystals 
[63, 64] could be enhanced in the same manner as the op- 
tical Hall effect in photonic crystals. It should be noted 
that the spin of the constituent particles is not always 
necessary; even a scalar wave can also have an internal 
rotation and a Berry curvature due to a periodic struc- 
ture breaking inversion and/or time-reversal symmetries. 



APPENDIX A: QUANTUM OR CLASSICAL? 

In this paper, we employ a quantum-mechanical for- 
malism in order to formulate our theory of a photonic 
system on an equal footing with that of an electronic 
system. The photonic system is described by an effective 
model where a dielectric medium is regarded as a classi- 
cal object. This quantization procedure of photons is cor- 
responding to that of electrons described by an effective 
model, e.g., a model with an effective mass-matrix and/or 
an effective (periodic) potential. In an electronic system, 
an electron is treated as a quantum object, and the equa- 
tions of motion for a semiclassical wavepacket are derived 
from an effective Lagrangian. In a photonic system, the 
counterpart of this effective Lagrangian is most naturally 
represented in a second-quantized formalism with keep- 
ing its close connection to that in an electronic system. 
This is because, in the Maxwell theory, wc; c;annot define a 
positive-definite probability density, while we can define 
a positive-definite energy density. (The second quantiza- 
tion is adopted to define the quantum wavefunction of a 
photon in Sec. 11 C, which cannot be directly represented 
by the field strength of an electromagnetic field.) How- 
ever, as long as we consider an approximately coherent 
wavepacket, the center of wavepacket coincide with the 
center of gravity as shown in Sec. 11 C. This fact enables 
us to link our quantum-mechanical theory and classical 
electrodynamics . 

One may wonder whether the optical Hall effect is 
quantum or classical. We cannon answer this question 
in a single sentence, but give some remarks on it as fol- 
lows. This effect comes from the particle-wave duality of 
an optical wavepacket and the geometrical/topological 
property of a wavefunction. Therefore, similar geomet- 
rical/topological effects are expected in various kinds 
of quantum/classical and microscopic/macroscopic wave 
dynamics, when wavepackets under consideration arc ap- 
proximately coherent. Indeed, we can formulate a theory 
for this class of phenomena based on a classical wave dy- 
namics of a macroscopic system, while its direct connec- 
tion to electronic systems is not necessarily clear. The 
confusion represented by the above question is mainly 
due to the situation in which we sometimes refer to 
wave equations for photonic systems as classical Maxwell 
equations and those for electronic systems as quantum 
Schrodinger equations, while both photons and electrons 
are quantum objects. The origin of this situation comes 



from two of advantages of photonic systems which lead 
to the success of classical Maxwell theory. In contrast to 
electronic systems, the statistics of photons is bosonic, 
and effective self-interactions between photons are usu- 
ally very weak. However, as long as we treat a photon 
and an electron in a single particle approximation, we 
can formulate both theories on an equal footing. 

It is beyond the scope of this paper to fix the termi- 
nologies "quantum" and "classical" common in photonic 
and electronic systems. Although we treat a photon as 
a quantum object through to the end of this paper, we 
refer to results obtained purely by wave dynamics of light 
as those obtained by classical electrodynamics. As long 
as we consider an approximately coherent wavepacket 
in a single particle approximation of quantum theory of 
photon or in a linear approximation of classical electro- 
dynamics, results obtained by both formalisms coincide 
with each other as shown in Sec. Ill A. 



APPENDIX B: ORTHONORMALITY OF 
EIGENFUNCTIONS 

The orthonormality, Eqs. (33a) and (33b) is approved 
with the orthonomality of Bloch functions, Eqs. (23a) 
and (23b). For example, Eq. (33a) can be shown by using 
the following relation, 

X I dre-^('=-'=')-C/„tfe(r)e(r)C/„^,,fc,(r) 

'^VEnkEn'k' 

/ dre-^('=-'=')-("+'-)[/„5,fc(r)e(r)C/„^,,fe,(r) 

^ Er ^k+En'k' ^i^E„,-E^,^,)t - k' + G) 

'^vEnkEn'k' ^ 

dr -^i''-'''yru^^^(r)r{r)U3,,Ar), (Bl) 



ws ^^ws 

where a represents an arbitrary lattice vector. Since the 
lattice momentum k and k' are in the first Brillouin zone, 
we can reach the result 



(*nAfe|e|^f'A'fe') 
Enk + En'k 



2\/ EnkEn'k 

Kn'hx'Hk - k'), 



e'^^-^-^^-nik - k'){U^^^\7\U^,y^) 

(B2) 



where we have used Eq. (23a). In the same manner, 
Eq. (33b), can be also proved by using Eq. (23b). 
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Next wc prove the orthogonality, Eq. (33c). From the 
definition of the inner product, Eq. (32), we can show 

= {En'k' - Er,k) j dr*f^fe(r,i)e(r)*f,;,,fe,(r,t), 

(B3a) 

(B3b) 

In the case of Enk = En'k', it is clear that Eq. (33c) 
is approved. Thus, in what follows, we consider 
the case of Enk 7^ En'k'- In this case, wc can 
easily show ^ dr^^^{r,t)*e{r)^^,yk'{r,t) = and 
/ rfr*^;^j,(r,t)/i(r)*^;^,^,,(r,t) = from the relations, 



iE^,k'-E^k) J dr^^^k{r,t)e{r)^^,yk'ir,t) = 0, 

(B4a) 

[Elk' -Elk) j dr^Sxk{r,t)12{r)Kx'k'{r,t) = 0. 

(B4b) 

Consequently, the orthogonality, Eq. (33c) is approved in 
all cases. 

The above relations are derived from the eigen equa- 
tions, Eqs. (22a) and (22b). For example, the relation 
for *^;^fc(T■,^) is proved as. 



El,k' J dr^^^kir,t)^{r)^^'yk'ir,t) 

= Elk j dr^^^kir, *) e (r)*^,,,^, (r, i), (B5) 

where *fl'^{r) = ^(r) is used in the transformation from 
the second line to the third line, and *e^{r) = e'(r) is 
used in the transformation from the third line to the 
fourth line. A similar relation can be derived also for 
^Sxkir,t). 



APPENDIX C: WAVEPACKET IN A PERIODIC 
SYSTEM 

Here we present details about an optical wavepacket in 
a periodic system. Basic features of the wavepacket arc 
discussed in Appendix CI. These features are helpful 
to understand the effect of an additional modulation su- 
perimposed onto a periodic structure, which is discussed 
in Sec. II C. Some comments on a gauge transformation 
in momentum space is given in Appendix C2. In Ap- 
pendix C3, we present detailed procedures to evaluate 
expectation values which appear in Appendix C 1. 



1. wavepacket 



We begin with the wavepacket defined by 



\W) 



BZ 



dkw{k,kc,t)^Zcxalxk\^), (Cla) 



w{k, fee, t) = Wr{k - kc)e-'^^*''^\ 



(Clb) 



where Wr{k—kc) is a real function, and Wr{k—kc) and Zcx 
satisfy the normalization conditions, dkw^{k — kc) = 
1 and |-2cAp = 1, respectively. We assume Wr{k — kc) 
has a sharp peak aroimd fee = j^idk uj^^k — kc)k. It 
should be noted that, rigorously speaking, we need to re- 
place this single photon wavepacket with a coherent (or 
squeezed) state wavepacket when wc apply the present 
formalism to a light beam with macroscopic number of 
photons. However, from the linearity of the Maxwell 
equations, the equations of motion for the single photon 
wavepacket is applicable also to the macroscopic coherent 
beam. 

In a fermionic system, we can define the position oper- 
ator as the center of the probability density of a fermion, 
which is positive-definite both in non-relativistic and rel- 
ativistic cases. However, for a relativistic boson, the def- 
inition of its position is nontrivial. In order to find an 
appropriate definition for the position of wavepacket, we 
firstly consider the energy and the position weighted by 
the energy density evaluated as follows, 

{W\H\W) ^ Er^k^, (C2) 
{W\n\W) ^ Er,kA^kJ{kc,t)-{z,\AnkJz,)]{C3) 

It should be noted that = in Eqs. (C2) and (C3) means 
that the above expectation values are evaluated under the 
assumption that the shape of w^{k — kc) is sufficiently 
sharp compared to the variations of Enk and A„fe(fe) 
around fee, and we neglected terms which depend on the 
shape of iv1{k — k^)- 

From Eqs. (C2) and (C3), the center of gravity is esti- 
mated as 

- ^kMkc,t)~{z,\KnkAzc). (C4) 

Comparing this result with the naive definition for the 
position of wavepacket, j0^wl{k — kc)'^k'&ik, t), we 
can reach the appropriate definition for the position of 
wavepacket, 

rc = dk wl{k - fee) [Vfe'!?(fe, t) - {Zc\A.nk\Zc)] ■ 

Jbz 

(C5) 

In order to check the property of the wavepacket, we 
consider the expectation values of physical observables. 
As shown in Appendix C3, the energy current and the 
rotation of energy current are evaluated as 

{W\V\W) ^ Enk^kEnk^, (C6a) 

{W\J\W) ^ rcy<Enkyk,Enk, + {Zc\Snk^Zc). 

(C6b) 
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where 



*nk 



XX' 



-'^[{-^kU^Xklxml^ 



^n^^kU^X'k) + {U^XklS X e -iS|C/„«v/c)] • 
I 



(C7a) 
(C7b) 
(C7c) 



It is noted that the first 
as the orbital rotational 
center of gravity, and the 
i.e. the rotation around 
for the locally isotropic 
arc scalar variables, e(r) 

the second terms in «S^j. and S"/^ are rewritten by using 
S X S = iS as 



term of (M^lJT'lVF) is interpreted 
motion, i.e. the rotation of the 

; second term as the internal one, 
the center of gravity. Especially 
system in which *e{r) and /u(r) 
and /i(r) , the contribution from 

H 



E 
Xk 



s X 7l-'s\u,fy^) 



(C8) 

This suggest that the internal rotation correctly includes 
the spin of the constituent particle, i.e., the polarization 
of light in the present case. However, it should be noted 
that the second terms of right-hand side of Eqs. (C7b) 
and (C7c) are not the whole contributions of spin. Ac- 
tually, when we consider the circularly polarized light in 
isotropic homogeneous media, all terms of the internal 
rotation give the same contribution and totally represent 
the rotation originated by the polarization. In addition, 
the internal rotation defined above contains the internal 
orbital one and the spin one generally. 



2. Gauge transformation 

When a system has a symmetry represented by the 
unitary matrix [Mnk]xx' , Maxwell equations are invariant 
under the transformation, 



\^nXk) 



]A'A|C^^A'fe)' 



(C9) 



where F = E or H. (Here we consider the case in 
which there are degeneracies indexed by the subscript 
A or A'.) For the sake of convenience, we call this trans- 
formation as the gauge transformation in fc-space. By 
this gauge transformation, A„fc, fink, ^nk and Snk are 



transformed as 

A„fe = M-^AnkMnk-iM-'^VkMr.k, (ClOa) 

^nk = M-j^UnkMnk, (ClOb) 

A„fe = M-jA„feM„fe, (ClOc) 

Snk = M-^SnkM^k- (ClOd) 

The gauge transformation of Bloch functions in Eq. (C9) 
is equivalent to that of the corresponding creation oper- 
ators as 



(Cll) 



In terms of this transformed operators, the wavepacket 
in Eq. (Cla) is represented by 



\W)= f dkw{k,K,t)Y,[M-l]x'xZcxaix'k\^). 

•'^'^ X,X' 

(C12) 



It should be noted that we have changed only the rep- 
resentation but not the physical state of wavepacket. 
Therefore, the expectation values of physical observ- 
ables, e.g., H, TZ, V and JT, must be gauge invari- 
ant. Indeed, we can easily show that the evaluations 
of H and P in Appendix C 1 are gauge invariant be- 
cause of the invariance of E^k- From Eqs. (ClOb)- 
(ClOd) and Eq. (C12), we can also show the invariance 
of {zc\Vtnk\zc), (2c|A„fckc), and {zc\Snk\zc)- However, 
it is not clear whether the evaluations of TL and given 
in Appendix C 1 are also the case. In order to confirm 
this point, it is enough to check whether the position of 
wavepacket Tc in Eq. (C5) is gauge invariant or not. In 
the representation of Eq. (C12), the derivative of phase 
factor Vfci?(A;, t) and {zc\Ank\zc) in Appendix C 1 are re- 
placed as 



Vk'^{Kt)^^{zcW''^''■''^Vk^ 



i-d(k,t) 



Zc) 



■i{Zc\ 



(fc,f)M 



nk 



nk 



(C13a) 

{Zc\Ank\Zc) {Zc\Mr,kAnkM-^\Zc). (C13b) 
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The above formulae and Eq. (ClOa) prove the gauge in- 
variance of as foUows, 



Jbz 



i{,Zc\ 



-{Zc\MnkKkM-l\Ze) 

rc+ dk wl{k - he) 
Jb 



IBZ 



xi{zc\ [M„fe(VfeM-fei) + (VfeM„fe)M-J] \zc) 
= r,. (C14) 

Combining this result and the gauge invariance of Enk 
and {zc\Snk\zc), we can confirm that the evaluations of 
and JT" in Appendix C 1 are also gauge invariant. 

3. Expectation values 

Here we presents the detailed evaluations of the ex- 
pectation values, i.e., the Hamiltonian H, the position 
weighted by the energy density 7?., the energy current V 
and the rotation of energy current JT, with respect to the 
optical wavepacket \ W) in a periodic system. The expec- 
tation value of an operator O is obtained from commu- 
tation relations between O and the creation/annihilation 
operators, 

= / dkdk'w*{k,kc,t)w{k',kc,t) 
Jbz 

dk dk' w*{k, kc, t)w{k' , kc, t) 



BZ 



x(0| 



|0), (C15) 



where we have introduced the abbreviation. 



nzck 



(C16) 



and this will be used also for the Bloch functions as 

lC5> = Jl'oA^nik)- (C17) 

The basic commutation relation between B{r) and D{r) 
can be represented in the following integral form, 

j dr dr' [*J (r) • B{r), D{r') ■ *2(r')] 

= -i j dr [Vr X *J(r)] • *2(r) 

= -i j dr ^l{r) ■ [Vr x *2(r)] . (C18) 



In particular, in a periodic system, by the above commu- 
tation relation and Eqs. (21a) and (21b), we can easily 
show that 



0"n\ki 



EnkEn'k' I dr 



^Txk{r)'^{r)^%X'k'{r) 

+*f;,(r)/J(r)*^,,,,,(r)^ 
= Enk5nn'5\w5{k - k'). 



(C19) 



In the transformation to the last line, we have used 

j dre-'^^-^'y-Fkk'ir) 

= V / dr e-*('=-'=')-("+'-)Ffefe. (a + r) 

= / — V5(fe-fe' + G)e-^('^-'^')'-Ffcfc,(r), 
Jws vws ^ 



(C20) 



where a represents an arbitrary lattice vector, and 
Fkk'{r) is a periodic function, i.e., Fkk'{a + r) = Fkk'{r) 

Since k and k' are in the first Brillouin zone, we have 
also used the following relation implicitly. 



Y^6{k-k' + G)e-'^''-'''>'- 



fe,fe'eist BZ 



= 6{k-k'). 



(C21) 



Then we obtain the result 



{W\H\W) 



Jbz 



dk wl{k - kc)Enk = Enk^- 



(C22) 



In the similar manner, we obtain the following commuta- 
tion relation which is needed to estimate the expectation 
value of Tt, 



ttnXk 

= E, 



n' X'k 

dr r 



*f,\(r)e(r)*f,,,fc,(r) 

+*^4(r)M(r)*^,Vfe'(r-)" 

{Vk~Vk')Kk-k') 
[{U^xk\*^\u5x'k-) + {U^xk\Wn'X'k')] (C23) 



= EnkEn'k' 
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In the transformation to the last line, we have used the 
relation, 



we obtain the result. 



J drre-'^''-'''y^Fkk'{r) 



Fkk'{r) 



dr 



2X 



ws 

(Vfe - VfcO 5^ -k' + G)e-'('=-'=')-'' 



(C24) 



and Eq. (C21). This commutation relation leads to the 

result, 

{W\TZ\W) 



= ^ dkEnk w*{k,kc,rc,t)'Vkw{k,kc,rc,t) 
^ Jbz l 

- [Vfeu;*(fe, fee, Vc, t)]w{k, fee, Tc, t) 

+ ^ / dkw'^ik- kc)E.ak 

= / dfeM;^(fc-fe,)£;„fc[Vfci9(fc,<)-(zc|Afc|zc)] 
Jbz 



Enk, [VkM^ct) - iZc\AkJZc)] ■ 



(C25) 



The expectation value of P is estimated by using the 
commutation relation. 



^nXk^ 



^1 '^il'A'fc' 

= EnkEn'k' j dr 



^n*Xk{r)^Kx'k'{r) 

= ^\/EnkEn'k S{k - k') 

X [{U^^k\iS\U^^x'k) - {U^xk\iS\U^'yk}] , 

(C26) 

where Eqs. (C20) and (C21) were used. Combining this 
commutation relation and the formula, 

Enk [{UgjiSlUg^,,) - {UgjiS\U,l^k)] 
= {U^.J [STl-'Pk -S + Pk- sj2-'s] \U^,^k) 
= {Ugj [v.sf] \U^,^k) 

= "^A, (C27) 



{W\V\W) 

= - / dkwl{k- kc) 
^ Jbz 

xEnk [(C/„lfeN^|f/lfe) - {uZjiSP^.^^}] 

= \ j dkwl{k-K)VkElk 
^ Jbz 

— Enkc'^kEnkc- (C28) 



The expectation value of is derived from the com- 
mutation relation. 



a-nXk 
= EnkE, 



-VEnkEn'k' [(Vfc - Vk')Sik - k') 
{U^Xk\iS\US^,^,) - {U^xk\iS\U^'X'k') 



(C29) 



where we have used Eq. (C24) in the transformation to 
the last line. It should be noted that, as in the previ- 
ous commutation relations, the above commutation re- 
lation is also restricted to the case in which both of fe 
and k' arc in the first Brillouin zone. In addition, our 
wavepacket is constructed from degenerate eigen modes, 
i.e, eigen modes with the same band index n. Thus the 
following formula, which will be proved later, is useful to 
estimate the expectation value of with respect to the 
wavepacket. 



Enk 



= -{Zc\^^nk\Zc) X {Enk^kEnk) + {Zc\Snk\Zc)- (C30) 



Combining the above relation and the commutation re- 
lation, we obtain the result, 
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{W\J\W) 



-I ■ 

Jbz 



dk w1{k — fee) 



[Vfei?(fe,i)] X (£„feVfe£;„fe) + ^ [{Ul^k\S X |VfeC/„«,^fe) 



BZ 



dfc W;^(fe - fee) [Vfet?(fe,t) - (^e|A„fe|^;e)] X {Enk^kEnk) + {Zc\Snk\Zc) 



rc X {Enk^kEnk) + (^^cl^nfecl^c)- 



(C31) 



The proof of Eq. (C30) needs basic but tedious cal- 
culations. Here we comment that the formula is con- 
firmed by using Eqs. (29a)-(31b), and the fe-derivatives 



of Eqs. (31a) and (31b). The outline of the derivation is 
given as follows, 



{Zc\A.nk\Zc} X {EnkV kEnk 
Enk_ \/ttE 

A \ . 



+ 



^^,JS X \VkU^,^^} + {VkUl^^l X S\U^,^^) - {U^^JS X \VkUZ^^) - (VfeCz.fel X s\uZ,^) 
(Vfet^^fcl X (VfeSf)|[/„\fc) + (Vfef/„«,^fe| X (VkE^WZ,^) 
(;7„\fc|(VfcSf ) X \VkU^,^^) + (C/,lfc|(VfcSf ) X |VfeC/,f,^;,)_ 

(Vfef/„lfe| X (e£4 - )| VfeC/„\fc) + (VfeCz.fel X OlElk " ^fe )|Vfe[/„lfe) 
'{UZMPk ■ S-e-'S X \VkUZ^k) + (VfeC/„lfe| X Pk ■ sT's\u^,^k) 
HU^Pk ■ SH-'S X |VfeC/„\;,) + {VkUZj X Pk ■ S*e-'S\UZ^k) 
I [(VfcC/„\fc| X {7El^ - Sf )| VfcC/„^,^fc) + {VkU^^J X (Mi?^fc - S^)|Vfct/„^,^;,) 

^ [(VfeC/„lfe| X (Pfe . Sri-'S - Sjl-'Pk ■ S)\U^,^^) - {U^^JiPk ■ Srt-'S - Srt-'Pk ■ S) X I Vfe[/^,^^) 
+(Vfe?7„lfe| X {Pk ■ ST'S - Sri-'Pk ■ S)\UZ^^) - {UZjiPk ■ Sri-'S - STfl-'Pk ■ S) x \VkUZ^k) 

= {zc\Snk\ze) + l^k X [{U^JiPk ' ^M^'^ " SjJ-'Pk ■ S)\uZ^^) + {U^JiPk ' Sl'^S - S^'^Pk ■ S)\U^,j' 

= {Zc\Snk\Zc). (C32) 



+- 



APPENDIX D: EXPECTATION VALUES IN A 
MODULATED SYSTEM 

When a modulation is introduced into a periodic sys- 
tem, the argument given in Appendix C are modified. 
Here we consider the modulation represented by Eq. (38) 
which is sufficiently weak and smooth. It is noted that 
the commutation relation Eq. (C18) are not modified 
even under any modulation. From this commutation 



relation, the creation and annihilation operators of ap- 
proximated eigen modes satisfy the same commutation 
relation as that in a periodic system as shown below. 
However, the approximated cigcn modes depend on the 
variable rc- Thus wc miist additionally take into account 
the operator '^ra'^nxk-r ^'-'^ ^^'^ derivation of the effective 
Lagrangian. Fortunately, we can show that the contribu- 
tion from ^rc^'liXk r vanishes by the following relation. 
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(Dla) 



1 
2 
1 
4 

= 0. 



7€('^c) 



V^^ In 
Vr„ In 



7/.(»'c) 



Enk — En'k 7 



VEnkEn'k 



(5(fc-fe')(f/nAfel^|f/„^A'fe) 



(Dlb) 



where, in each commutation relation, we have used 
Eqs. (21a) from the first to the second expression, and 
Eqs. (C20) and (C21) from the second to the third ex- 
pression. 

Next the expectation values of H and TZ in modulated 
system will be estimated. In the effective Lagrangian, 
we retain up to the first order with respect to the deriva- 
tive of 7e(r) or 7/i(r). Therefore we need to estimate 



the expectation value of H up to the first derivatives of 
the modulation functions. As for the expectation value 
of IZ., we may neglect derivative terms as was discussed 
in Sec. II C. The commutation relations needed to esti- 
mate {W\H\W) and (WITZIW) are calculated by deriva- 
tive expansion with respect to the modulation functions. 
Firstly, the commutation relation for (VF|i/|VF) is esti- 
mated as follows. 



--E, 



nk;rcEn'k';rc J dv 



^2(^ .^n*xkirmr)^^'yk'ir) + z^^^^xkirMr)Kx'k'ir) 



1 



— 2 ^nk;rc^n' k;r, 



^^g-i(fe-fe')-r-= 
+ ••• . 



~6{k - fc')[(C/^Afc|e|C^,?A'fc> + (U^xklW^'X'k) 
(Vfe - Vk')~S{k - k')] ■ \[Vr^ ln7e(re)] {U^xkl'^lU^'X'k') + [Vr. In7^(re)] {U^xkmU^'X'k') 



(D2) 



The first and second terms in the second expression come 
from the terms of zero-th and first order of (r — Tc) re- 
spectively. In the transformation to the last expression, 
we have used Eqs. (C20), (C21) and (C24). 



In the same manner, restricting to the case in which 
both of k and k' are in the first Brillouin zone and using 
Eq. (C24), the commutation relation for (W^|72.|W^) is 
estimated as follows, 
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= E„k 



Tc-^n'k' \rc j '^f ^ 



^nxkir)- e(r)*f,;,,fc,(r) 



^nxkir)- nrX'X'k'ir) 



+ • 



(D3) 



I 

From these commutation relations, we obtain the fol- of gravity as (VF|7?.|M^)/(W|iJ|H^) = Vk^{kc,rc, Zc,t) ■ 
lowing results which leads to the estimation of the center {zd^nkjzc), 



{W\H\W) 



{W\-R\W) 



dk wl{k - kc)Enk;r 



BZ 



Vrcln[7e(rc)7^(rc)] • V k'd{k,rc, Zc,t) - 



[Vr., In7e(rc)] • {zc\-f^^kW) + [V^e ln7/.(rc)] • (^c|A^fe|^c) + • 

In [le{rc)ln{rc)\ ■ (2c|A„fcJzc) 
[V^, ln7e(re)] • (^clA^fcJ^e) + [v., In7^(re)] • {z,\A^k^\z,) 



Enk 



Vr^ In 



{Zc\AnkJZc) 



E, 



BZ 

+ I dk Wi{k - K)Enk;r. I £ I VfeiJ^.^^) - (Vfc[/^,^fe| e |C/^,^fe) 



(D4a) 



dk Enk;rc W* (fe, fee, Tc, Zc, VfeW(fe, fee, Tc, t) - [VfcW*(fe, fee, Te, Z^, t)]w{k, fee, Te, t) 



JBZ 



+ ■■■ 



dkwl{k-kc)E,^k-ra ^k'&ik.rc.Zc.t) - {zc\^^nk\zc) 



= EnkcWo \Sk,'&{kc, Vc, Zc, t) - (2c| A„fc^ l^c)] . 



(D4b) 



r 



The above estimation for the center of gravity suggests 
that the position Vc defined by Eq. (42) may be regarded 
as the center of gravity even in the case with modulation. 
In the derivation of the effective Lagrangian, we need to 



estimated the inner product between the wavepackct and 
its time-derivative. Finally we present the detail for the 
calculation of this product, by regarding Eq. (42) as the 
definition of the center of wavepacket. 
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= i [ dkw*{k,kc,rc,Zc,t)-^w{k,kc,rc,Zc,t) + i{zc\zc) 
at Jbz 

+i / dk dk'w* {k,kc,rc, Zc,t)w{k' ,kc,rc, Zc,t){0\a„^^k-,- 



''c nzck';rc 



dkw^{k - kc)—i9{k,rc,Zc,t) +i{zc\zc) 



-j 

Jbz 

= -kc- I dk[Vk^w'^{k-kc)]i!}{k,rc,Zc,t)+i{zc\zc) + -^ [ dkwl{k - kc)^{k,rc,Zc,t) 
Jbz Jbz 

= kc - [ dk[VkW^{k-kc)]'&ik,rc,Zc,t)+i{zc\zc) + [ dkw^{k - kc)'&{k,rc,Zc,t) 
Jbz dt Jbz 



— kc ' 

— kc • 



BZ 



dk wl{k - fee) [Vfcl?(fe, Tc, Zc, t)] + i{Zc\Zc) + ^ / c''^ "'rC^ - fec)^(fc, ■Tc, -Zc, 



BZ 



rc+ dkw^{k - kc){zc\Ank,\z, 
Jbz 

fee • re - fee • (^e|A„fej2;e) +i(2e|ie) + ^ 

r 



+ «(-2e|4) 



dk 'w1{k — fec)'!?(fe, Tc, Zc, t) 



Jbz 



BZ 



dk wl{k — kc)^{k,rc, Zc,t) — kc - rc 



(D5) 



where Eq. (Dlb) and anz^k;rjO) = are used in the fe-derivative in the above expression in terms of the 
transformation from the first expression to the second Feynman-Hellman relation. However, in systems with 
expression. the gauge symmetry, even if \U^-^^) is a Bloch function 

of a physical state, its derivative may have an unphysical 
component proportional to \K). In other words, Bloch 
functions and their derivatives should be expanded by 
non-orthogonal bases as follows, 



APPENDIX E: BERRY CURVATURE AND 
INTERNAL ROTATION 

In a system with generic periodic structure, it is a 
tough work to analytically calculate the Berry curvature 
and the internal rotation. However, it is easy to obtain 
them numerically by rewriting inner products of Bloch 
functions and their momentum-derivatives to the prod- 
ucts of conventional expectation values. Here we present 
some formulae which are convenient for numerical calcu- 
lations. 

For latter convenience, we separate the Berry curva- 
ture as, 

0„fe = i [f2^fe + n^fe] - iA„fe X A„fe, (Ela) 



+ E mT^]-\G,G'){K'\7\Vk}{E2^) 

G,G' 



E \K)[T^]-\G,G'){K'\Jl\Vk){E2h) 



G,G> 



where rf(G,G') = K*eK' and r^(G,G') = K*jlK'. 
(Elb) By using the above expansion and the Feynman-Hellman 
relation derived from Eqs. (31a) and (31b), we can 
where F = E ov H. In the following, we rewrite the rewrite the Berry curvature as 



fo-E 1 _ (^TtAfcl [^fc"fc] l^mA"fc) ^ (^mA"fcl [^fc'^fc] l^nA'fc) , /rrE \*^!rE]-l q^\ttE \ (-p'i<.\ 

\^t^i^ _ ■ ^n\i^ [^fc-fe] l^mA"fc) ^ (^mA"fcl [^fc-fe] \^n\'kl , irjH Ittrr^l -1 et?! rr-H \ /pou^ 

{^nk\xx' = [{U^xk\S\U^x'k) + {U^xk\S\U^X'k)] ■ (E3c) 
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Thus {i = E, H) is enhanced when the band comes close to other bands in energy, with the enhancement being 
inversely proportional to the square of energy difference. In contrast, A„fc does not have such an enhancement. 
Though Eq. (E3c) seems to diverge at Enk — » (fe — > 0), it is not the case, as we will see in Sec. IIIC for a specific 
case. In the long wavelength limit A; ^ 0, the propagating light becomes insensitive to spatial modulations of yu(r) 
and e(r), and the medium is regarded as uniform. Because A„fc = for a uniform isotropic media, a generic periodic 
medium in a long-wavelength limit also shows A„/e 0. 

In the same manner, the internal rotation is also rewritten as follows, 



XX' 



XX' 



.X"kl 



E I 
mX"k\ 



\ttE \ 
^nX'k/ 



E^^ -E^^ 
nk mk 



+E\ 



{U^xu\ 



[rf] 



X'kl 



{U^xk\ [Vfe^f] \Kx"k) X {U^^x"k\ [Vfe^f] W^x'k) 



T?2 

mk 



+El,,{U^^k\Wkr^Sll\UnX'k) - i{U^xk\S X e-iS|C/f,,fc) 



(E4a) 



(E4b) 



It should be noted that ^^k ^nk li^ve very similar 
expressions to <S„fc and <S„fc, respectively. This suggests 
that there are always some kind of rotation when the 
Berry curvatures are nonzero. In this sense, we have 
generalized the argument for the quantum Hall system 
in Ref. [34] to a photonic system. In the quantum Hall 
system, the internal rotation is the internal orbital rota- 
tion originated by the cyclotron motion under an external 
magnetic field. On the other hand, in the present case, 
the internal rotation are the combination of the polar- 
ization and the internal orbital rotation originated from 
a periodic structure. When there is no periodic struc- 
ture, anisotropy nor inhomogencity in ?(r) and 
Eqs. (E3a)-(E4b) are reduced to the Berry curvature, 
p-f73, and the spin divided by e/x, ^ • ^0-3. These con- 
tributions come only from the terms including the spin 
operator S, and A„fc — 0, i.e., nonzero A„fc is originated 
by the anisotropy or the periodic structure of e* and /I- 
Even in generic cases, A„fe has a unit of a length, and 
its magnitude is a lattice constant at most. 



APPENDIX F: TRANSVERSE SHIFT IN 
CLASSICAL ELECTRODYNAMICS 

Here we prove the consistency between our result for 
the transverse shift (Eq. (63)), which is consistent with 
the TAM conservation for individual photons (Eq. (65)), 
and the result by Fedoseev [26, 27], which is based 
on classical electrodynamics. In Refs. [26, 27], each 
wavepackct is constructed as a superposition of plane 
waves with wavevectors k = kc + k, where k distributed 
around zero vector. (In the notation of Refs. [26, 27], fee 
is represented by K.) The polarization vector of each 
constituent plane wave is defined by Eq. (23) in Ref. [26] 



with Eq. (7) in Ref. [27], 

e«(K) = zi'^iK)s'^^\K) + zl'\K)p<-^\K), (Fla) 

where j — i,p,T for incident, reflected and transmit- 
ted beams, respectively, 2;P^(k) and Zp\k) (|2;P^(k)|^ -|- 

\zjP {k)\'^ = 1) represent the polarization state of each 
plane wave, s'^^\k) and p'^^\k) are the s- and p- 
polarization vectors deflned by 

^"'W^^Wxi, (F2) 

where n = (0,0,1) is normal to the interface, and we 
consider the same configuration of the interface and 
beams as those in Sec. Ill A. The relation between 
the present notation and that in Ref. [27] is repre- 
sented as zi^\K) ^ AW(/«)/V|A0H«)P + |SW(k)|2 
and zI^\k) ^ B(^)(k)/V|^W)(k)|2 + |BW)(k)|2, n ^ 
N and fc/|fe| ^ m{K). 

By the Maxwell equations, ^s'''^^(k) and Zp'''^\k) are 
exactly given by 

^iH-) = I , (F3a) 

^|il^)(K)^W|2-H4-')(K)4^)|2 

4^')(k) = t^^\K)z^\K) ^p^^^ 

V|ii^'^('c)4'^|2+|4^)(K)4*)|2 

where j = p or r, t'f^^K) and t^\K) are the ampli- 
tude reflection coefficients for the s- and p-polarizcd plane 
waves, ii^ ■*(/«) and t^\K) and the amplitude transmis- 
sion coefficients for the s- and p-polarized plane waves, 
i.e., fSt^ ^ tlf^ ^ Rp, t^s^ ^ Ts and 4'^^ ^ Tp in our 
notation in Sec. HI A. 
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In our constitution method for an incident wavepacket, 
the polarization state of each constituent plane wave, 

i.e, the set of z'i^ and Zp \ is independent of k,. Oth- 
erwise, the concept "an clliptically-polarized incident 
wavepacket" gets fuzzy (see Sec. Ill A 3). Thus, this is 
a natural definition for an elliptically-polarized incident 
wavepacket. Its polarization vector is represented also in 
the following form. 



(F4) 



where m is a complex constant, representing the polariza- 
tion state. This m is identical with m defined by Bliokh 
et al. [50], and related with our \z^) in Sec. Ill A by 



1 — im 



It yields 

iz'\T\z') = — [1 - Imp, 2JR(m), 23(m) 



(F5) 



(F6) 



which is used for comparison between the results here 
and those based on our theory of the TAM conservation 
for individual photons. 

Wc now calculate the transverse shift from Eqs. (15)- 
(17) in Ref. [27]. The result is a sum of two terms 



dy 



(F7) 



where j = p or r, and 5y^''^ <-> 6y^ and Sy^"^^ Sy^ 
in our notation in Sec. Ill A. From Eqs. (13a), (13b), 
(17) and (18) in Ref. [27], the second term of right- 
hand side, h^^^\ is proportional to the K^-derivative of 

3[lnzp-'(K) — In Zp\k)] at /« = 0. The amplitude re- 
flection/refraction coefficients depend only on the polar 
angle, and thus their derivatives by Ky at k = are zero, 
because the y-component of kc is zero in the present con- 
figuration. As was mentioned previously, zi^^ and Zp^ are 
independent of k. Therefore, from Eqs. (F3a) and (F3b), 
h^"^ = {j = p, r), and we have 



Sy 



+i 



,n- [e^^Ho) X e(^:'*(0)] 

|n X fcO)| 
n- [e«(0) X e«*(0)] 
|n X feW| ' 



(F8) 



where k^^^ {j = p or r) are mean wavevectors for re- 
flected (p) and transmitted (r) wavepackets. Note that, 
the correspondence between these wavevectors and those 
in Sec. Ill A are k^'^ ^ k^ , k^^^ ^ k^ and fe(^) ^ k^. 
Eq. (F8) is identical with Eq. (63), showing an equiva- 
lence between Fedoseev's theory based on classical elec- 
trodynamics and ours. 



Finally, wc rewrite Eq. (F8) in terms of our notation 
in Sec. Ill A. For partial reflection, Ap and ^l^ are real, 
and we get 



Sy' 



23(m) 
sin 6 1 



{As /Ap ) cos 9 a cos 9 1 
l + iAJApf\m\^ ~ l + \m\^ 



(F9) 



where A = T ot R. By rewriting Eq. (F9) in terms of 
\z^), the shift is equal to our result in Eq. (68) but not 

to Eq. (5) in Ref. [50]. For total reflection, Rp and R^ 
are complex numbers with = \Rs \ = 1, and wc get 



6y^ = 



-2 cos 9i 9(TO)[5i(i?;i?,) + 1] + 3?(m)9(i?;i?s) 



sin 9 1 



1 + |m|2 



(FIO) 

This is exactly the same as ours in Eq. (69). To sum- 
marize, for every case, the calculation based on classical 
electrodynamics gives the identical result with ours based 
on our quantum-mechanical formalism, and this result is 
consistent with the TAM conservation for individual pho- 
tons. 



APPENDIX G: BERRY CURVATURE IN A 
TWO-DIMENSIONAL PHOTONIC CRYSTAL 

In order to discuss the TM and TE modes, it is conve- 
nient to introduce the following unit vectors, 

K Cz X ex 

e/f = — , ei = 



\ez X Cif 



and the Bloch functions are represented by 



for the TE modes and 



^TM mfc/ 
M|C^TM mfc) 



Bz (S) IC^TIvI mfc)) 



(Gl) 



(G2a) 
(G2b) 



(G3a) 
(G3b) 



for the TM modes. The superscripts, D and B, mean 
that they correspond to the electric and magnetic flux 
densities, respectively, satisfying the transversality con- 
dition, i.e. being perpendicular to K. 

The matrices for the eigen equations given in Eqs. (26c) 
and (26d) are simplifled in the case with scalar e(r) and 
/u(r) as follows, 

sf(G,G') = e{G,G')^i-\G,G'), (G4a) 
S^(G,G') = e{G,G')e-HG,G') (G4b) 

where QkiG, G') = {K-K'I-K'®K), and their deriva- 
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tives are represented by because G and G' have no 2;-components. 



-Ky ^k. 




Vfc,efe(G,G') = I -K'y K,+K', 

-k, K,+K'^ 



\7k,ek{G,G') 



From the above formula together with Eqs.(E3a), 
(E3b) and by setting kz = 0, we can easily show that 
the Berry curvature of a non-degenerate TM (TE) mode 
has only ^-component, 



(-^TM rife -^TM mk 



(^H,z _ r, \ " ^ [(^TMwfcl [^fc:r^fc] I^TMmfc)(^TMmfcl [^fc»"fc] I^TMrafc)] , /rrH \ , trH]-! Qz , \ttH \ 

y-'^TM nk ^mTM k) 



(G6b) 



)E,z _ r, \ ^ ^ [(^TEwfcl [^fc»"fc] I^TEmfc)(^TEmfcl [^fc^'^fc] I^TEnfc)] 



r~,E,z _ o \ ^ ^ LN^TEwfcl ["fex^fcj I^TEmfc/X^TEmfcl Ky^k j I^TEnfc/J , /ttE \^\-rE]-l qz^\ttE \ /nRr.\ 

^^TEnfc - ^ r52 Zrp2 \2 1- \'^TEnfc|e[i fcj eFxEnfe/' l^OCj 

TO^„ l-^TE nfe -'^TE mfe-' 

_ o \" ^ [(^TOnfel [^fcx-fc^] \Ute mk) i^TE mk\ [^ky^k] Wte nk)] lnaA\ 

^'TEnfe - ^ rP^ Zrp2 \2 • [^M) 

V-^TEnfe ^TE mk) 



r 



Note that the Berry curvature does not necessarily de- 
crease with the energy increases in contrast to the sys- 
tem without periodic structure, because nearly degener- 
ate points due to the band structure enhance the magni- 



tude of the Berry curvature. 

Following the same argument as the Berry curvature, 
the internal rotation also has only the z-component for 
non-degenerate bands. 



qE,z 
^TM nk 



H,z 
TM nk 



cE,z 
'TE nk 



crH,z 
^TE nk 



m^n 



^TM nk 



^TM mk 



p2 

-^TM nk 



p2 

-^TM mk 



E 

m^n 

+ 2 [-^TMnfe(^TMnfelMPfe^] ^ f^-l^TM nk) i^TM nk\^ ^'^'^I^TM nfe)] 
^ [(^TEnfel 



[Vfc.Sg] \U^^r^k){U^Emk\ [Vfe.Sg] It/fEnfe)] 



m^n 



p2 _ p2 

-^TE nk ^TE mk 



+ 9 [-^TEnfe(f^TEnAfcl^[rfe] ^"S'^^I^TE nA'fc) + (^TE nfe IM ^"^^If^TEnfe)] ) 



E 



^ mEnk\ [V..S^] \U^^^k){U¥Emk\ [Vfe^Sf] \U^^ „k) 



p2 

-^TE nk 



p2 

-^TE mk 



(G7a) 



(G7b) 

(G7c) 
(G7d) 



It is interesting that the internal rotation of a photon can be perpendicular to the propagating direction in a 
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two-dimensional photonic crystal. 

We also calculate A„fc from Eq. (E3c). For non- 
degenerate bands, there is no contribution to the Berry 
curvature fink from the vector product of Axm(te) nfc i^i 
Eq. (Ela), because Axm(te) nit is a simple vector vari- 
able, not a set of matrices. Meanwhile, Axm(te) nk may 
modify the energy spectrum when a modulation is ap- 
plied. A„fc is given as follows 



1 



^TM nk 



^TE nk 



^-^TM nk 



^^TE nk 



(G8a) 

^[{U^Enk\t^-'Pk^-'f^-'\U^Bnk)]- 

(G8b) 



In many cases, we can approximately regard the magnetic 
permeability fj, to be constant. Then AxMnfe vanishes 
from Eq. (G8b), whereas AxEnfe does not in general. 



APPENDIX H: REMARKS ON Ate nk 



to a vector potential and a magnetic field.) In such a 

case, Ate nk and IIte nk are not necessarily minor cor- 
rections. These corrections may become enough measur- 
able for a generic modulation additional to a periodic 
structure, while, for assuring the validity of our argu- 
ment, we mainly consider a slowly-varying modulation 
in this paper. 



(a) TE 1st [A' /a, A' /a] logdAli) 




ka 



(b) TE 2nd [A'/a, A'/a] 





ka 



FIG. 6: Difference between the electric and magnetic parts of 
the Berry connection for each of (a) the TE first band and (b) 
the TE second band, i.e., Ate nfe ~ {-^i 



'■TE nk 



-^TE nk)/^- 



Here we evaluate Ate nk for the two-dimensional pho- 
tonic crystal discussed in Sec. Ill C and see its effect on 
the energy dispersion and group velocity for each of the 
first and second bands of TE mode. From Eq. (48), an 
additional correction appears in the energy of each TE 
mode as 



^TE nfcc;r-c 



E' 



TE nkc-.Va 



l-[V,Jn7,(r,)]-ATE„fc. (HI) 



where i?TEnfc,;r, = 7€{rc)ETEnk,- Figure 6 shows 
Ate nk for the first and second bands of TE mode and 
we can see Ate nfe ^ 0.1a. Therefore, the correction is 
at most a few percent as long as the modulation is suffi- 
ciently weak, i.e., jaV^^ \n-f^{rc)\ <C 1. In order to make 
the argument complete, we also calculate a correction to 
the group velocity of a TE mode, 

= ^ ka^TE nka;ra " ^fec [^r<,7e(''c)] ' AtE nfc^ -^TE rife^ 
— Vfc^i?TEnfee:r<, + HTEnfej^fec, (H2a) 
n^E nk = ^l^^E nk + [^l 1^ ^^TE nk] A^^E nfe- (H2b) 

Here we used the relation kc = — [Vr^7c(rc)]i?TE nfc^ for 
smooth and weak modulation. By plugging Eq. (H2a) 

to the equation of motion for Tc in Eq. (53a), HTEnfc 
is a variable to be compared with the Berry curvature. 
Figure 7 shows that the effect of HTEnfc is negligibly 
small compared to the effect of the Berry curvature in 
the present case. However, it is noted that, even when 
a Berry connection is nonzero, the corresponding Berry 
curvature can vanishes. (This is easily understood by 
the analogy of a Berry connection and a Berry curvature 



(a) TE 1st [n"/a , W'la ] iog(in'l/a-) 




(b)TE Istjn-'/a^ n"/a1 i„g(in'IA,'-) 




ka 



(d) TE 2nd [W'la, W'la] 



'g(in'l/a') 




ka 



FIG. 7: IIte nfe, which is related to the correction of group 
velocity for each of (a) the TE first band with a a;-directional 
modulation, (b) the TE first band with a y-directional mod- 
ulation, (c) the TE second band with a a;-directional modu- 
lation and (d) the TE second band with a y-directional mod- 
ulation. 



APPENDIX I: DIFFERENCE FROM 
MAGNETICALLY INDUCED DEFLECTIONS 

It was proposed theoretically [39] and observed exper- 
imentally [40] that, in a Faraday-active random medium 
subject to a magnetic field perpendicular to an incident 
beam, the diffusion flow of light is deflected in a direc- 
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tion perpendicular to both the incident hght beam and 
the externally applied magnetic field. This effect seems 
to be more similar to the conventional electrical Hall ef- 
fect than the optical Hall effect is, because the effect is 
caused by the external magnetic field and the direction 
of deflection is perpendicular to it. However, it should 
be noted that, unlike electrons, photons are not charged, 
and their orbital motions do not directly couple to an 
external magnetic field. This effect is theoretically inter- 
preted by the magnetically induced off-diagonal compo- 
nents of a diffusion tensor and experimentally proved to 
be due to the magnetically induced changes in the op- 
tical properties of scatterers [39, 40]. In this sense, this 
effect is similar to the anomalous Hall effect due to the 
skew scattering mechanism, rather than to the conven- 
tional Hall effect. On the other hand, the optical Hall 
effect is originated by the anomalous velocity of an opti- 
cal wavepacket which appears without external magnetic 
field nor scatterers. 

This kind of phenomena, i.e., magnetically induced de- 
flection, is not restricted to random media. The deflec- 
tion of light by a magnetic field in a nonscattering ho- 
mogeneous medium has also been discussed theoretically 
[41] and observed experimentally [42]. When the effect of 
absorption in a Faraday-active medium is negligible, the 
linear effect of external magnetic field B on this medium 
is described by the dielectric tensor. 



n^{Sij + 2i€ijkAk), 
^ R 



(Ila) 
(lib) 



where n is the refractive index of the medium in the case 
of B = 0, 3?7 and ^7 represent the strength of the mag- 
netic circular birefringence and that of magnetic circular 
dichroism respectively, while we set ^7 = 0. The eigen 
modes of the dielectric displacement D = *eE in such 
a medium are explicitly given in Ref. [43], and they are 
represented in terms of the orthogonal unit vectors e^, 
ee, and in the spherical coordinate of the fe-space as 

D+ oc (es • ek)eg + i{CB + A[e_B x efe[^)e0,(I2a) 
D_ oc {Cb + A[eB X ek\^)ee - i{eB ■ eA;)e^,(I2b) 



Cb = \/{&B ■ efe)2 + A2|eB x ek\ 



(I2c) 



Here es is a unit vector defined by A = A es with the 
condition cb • Cfe > 0, and || eg x ej,, e^i = x efe. 
(When 37 7^ 0; A is a complex- valued parameter.) These 
eigen modes have the dispersion relations and the group 
velocities. 



vk 



Vl-2A2|es xefe|2T2ACB' 



k 



Cfc T T^les X efc|(eB • ek)ee 
Ob 



(I3a) 
(I3b) 



where v = 1/n. The direction of Poynting vector of each 
mode coincide with Vk± as long as ^7 = 0. It should 



be noted that the deflection occurs within the plane de- 
termined by k and B. The angle 69 between the propa- 
gating directions of two eigen modes with the same k is 
given in Ref. [43] and represented in the present notation 
as 



66 = 2arctan 



A|eB X efc|(eB ■ e^) 
V'(eB • efe)2 + A2 [cb x e^l 



(14) 



For the exact Voigt geometry [es-ek = 0), there appears 
no deflection [43, 44]. The physic of this phenomenon is 
intuitively interpreted by considering the first order per- 
turbation with respect to A and the situation in which 
the angle between cb and Cfc are not close to the Voigt 
geometry, i.e., bb ■ e.k ^ |A||eB x e^l^. The approxi- 
mated eigenvalues and group velocities are represented 
as follows, 



V±,k 



v{k±A-k), 
v{ek ± A). 



(I5a) 
(I5b) 



This effect comes from the magnetically induced change 
in the dispersion relation of each mode due to the 
Pitaevskii magnetization, ±t;A • k [41]. On the other 
hand, the optical Hall effect is caused by the anomalous 
velocity due to the geometrical propriety of a wavepacket. 

In the above perturbative picture, D± are approxi- 
matcily ciquivalent to right/left circularly polarized modes 
which have the spin angular momenta, ±6^. Therefore, 
the above interpretation based on the Pitaevskii magne- 
tization means that an external magnetic field couples 
to the spin of photon through a Faraday-active medium. 
Prom this consideration, we reasonably expect that an 
external magnetic field couples not only to the spin but 
also to a generic internal rotation of photon in the form 
of dipole coupling. (Consequently, this effect is expected 
for Laguerre- Gauss beams [59] which have internal or- 
bital angular momenta.) As shown in Sec. IIIC, there 
appear eigen modes with large internal rotations in a two- 
dimensional photonic crystal without inversion symme- 
try. Here we take the configuration in which the photonic 
crystal is periodic in the .xy-planc and uniform along the 
z-direction. Considering an eigen mode with kz = 0, 
its internal rotation is oriented in the ^-direction, i.e., 
perpendicular to its propagating direction. Thus, when 
the photonic crystal is composed of Faraday-active me- 
dia and subject to an external magnetic field, it is ex- 
pected that the magnetically induced deflection can be 
enhanced. In addition, this effect would be observed even 
in the Voigt geometry. The details of this problem is 
beyond the scope of the present study and we will dis- 
cuss it elsewhere. Here we just note that this effect in 
a Faraday-active photonic crystal is due to the magnet- 
ically induced change of dispersion relation as well as 
that in a homogeneous Faraday-active medium, and is 
different from the optical Hall effect in a photonic crys- 
tal discussed in Sec. IIIC. 
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